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Abstract 


This users manual provides in-depth information concerning installation and 
execution of Laura, version 5. Laura is a structured, multi-block, compu- 
tational aerothermodynamic simulation code. Version 5 represents a major 
refactoring of the original Fortran 77 LAURA code toward a modular struc- 
ture afforded by Fortran 95. The refactoring improved usability and maintan- 
ability by eliminating the requirement for problem-dependent re-compilations, 
providing more intuitive distribution of functionality, and simplifying inter- 
faces required for multiphysics coupling. As a result, LAURA now shares 
gas-physics modules, MPI modules, and other low-level modules with the 
Fun3D unstructured-grid code. In addition to internal refactoring, several 
new features and capabilities have been added, e.g., a GNU-standard instal- 
lation process, parallel load balancing, automatic trajectory point sequencing, 
free-energy minimization, and coupled ablation and flowfield radiation. 
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1 Introduction 


The users manual consists of seven sections. Section 2 gives an overview of 
new features, capabilities, and bug fixes. System requirements and installation 
are covered in Section 3, followed by code execution instructions in Section 4. 
Section 5 presents input files, their formats, and detailed information on their 
contents while Section 6 covers output files. Ancillary utilities are explained in 
Section 7, and the last section, Section 8, presents illustrative example cases. 
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2 New in This Version 


Laura v5.2 offers several new enhancements and bug fixes since the previous 
released version, v5.1 [1]: 

• Major Enhancements 

— Automated radiation mechanism selection and streamlined user in- 
terface 

— Added multi-processor per block capability for coupled radiation 

• Minor Enhancements 

— Extended laura_blayer.dat to accommodate i and j surfaces 

— Added Troubleshooting and Support sections 

— To streamline the assignment of Plot3D function file variable 
names, Tecplot™ .nam files are now produced instead of macros 

— To facilitate plotting multiblock cases, a FieldView™ .fvbnd file is 
produced for volume data 

— To play nice with *nix, codes now supply non- zero exit code upon 
failure 

— To improve quality assurance, added technical debt build, unit test 
build with code coverage statistics, equilibrium air regression tests, 
and Shuttle Orbiter regression test 

— To facilitate configuration management and code support, laura 
now emits version string via -V or — version command line options 

• Bug Fixes 

— Stop gracefully when tdata not found 
— Multiprocessor run hangs for invalid namelist 
— Equilibrium air data files not installed 
— Intermittent failures of coupled ablation cases 
— Load balancing for large numbers of processors 
— Out of bound array in coupled ablation cases 
— Logic error in algebraic turbulence models 
— Pseudo-cell assignments across block boundaries 
— Unrelated blocks exchanging information 
— Out of bound array for wake cases 
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— Undefined memory usage in radiation 
— Out of bound array for radiation 
— Initialization of ablation surface properties 
— Full-fraction line implicit solves for certain multiblock grids 
— Grid alignment for cases with axis singularity 
— Partial-fraction line implicit for 3D cases 
— Non-pivoting option for linear solver 

For a complete revision history, please see the ChangeLog hie in the source 
distribution. 
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3 Installation 

Laura requires a Fortran 95 compiler, and if parallel processing is desired, 
a Message Passing Interface (MPI) implementation. 1 Some optional utilities 
require Ruby. 2 The installation and subsequent execution of Laura assumes a 
Unix-like operating system or compatibility layer. 3 After the code is unpacked 
from the LAURA release tarball, 

7, tar zxf laura-5.2-Z.tar.gz (unpack gzipped tarball) 

7» cd laura-5 . 2-Z 

where Z is a revision track number. Laura is installed via GNU build sys- 
tem, 4 which entails executing a sequence of four commands: configure, make, 
make check, 5 and make install. 6 

3.1 Sequential installation 

To configure, compile, test, and install a sequential version of Laura for use 
with a single processor, first make a subdirectory of laura-5. 2-Z to store the 
configuration. For example, 

7o mkdir g95-seq 
7» cd g95-seq 

if using the g95 Fortran compiler; 7 and then proceed with the typical GNU 
build sequence, 

7» ../configure FC=g95 — prefix=$PWD 
1 make 

7, make check 
7. make install 

Note that configured --prefix option specifies the root directory for in- 
stalling build artifacts — the default is /usr/local. In this example, it is set 
to the current working directory, g95-seq so executables will be installed 
in g95-seq/bin and data files will be copied to g95-seq/share/laura and 
g95-seq/ share/physics_modules directories. 

To use Laura and associated utilities, set your search path to include 
$PWD/bin, e.g., 

setenv PATH ${PWD}/bin : $PATH (for csh) 
export PATH=${PWD}/bin : $PATH (for sh) 

1 For example, OpenMPI or MPICH. 

2 See ruby-lang.org. 

3 For non-Unix-like systems, compatibility layers are available from mingw.org (minimal) 
and cygwin.com (maximal). 

4 See gnu.org/ software / autoconf/. 

5 The make check command is optional. It will attempt to run small test cases. 

6 The make install command may require administrator privileges depending on your 
installation location. 

'g95.org 
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3.2 MPI Installation 


An MPI-enabled installation (to allow multiple processors) is similar to the 
sequential installation except the configuration command has the --with-mpi 
option instead of the FC Fortran compiler variable, e.g., 

7. ../configure — prefix=$PWD \ 

— with-mpi=/usr/local/pkgs/ompi_l . 2 . 8-intel_ll . 0-028 

Another difference is that an MPI-enable configuration will produce an exe- 
cutable named laura_mpi instead of laura. 

In either case, config.log contains a record of the configuration com- 
mand used, and conf igure’s --help option details all available configuration 
options. 
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4 Execution 

The following steps outline a typical simulation cycle. 8 

Step 1. To start Laura, a Plot3D structured grid file is needed — see 
Section 5.2 on page 14 for more info. You may externally generate a 
grid using grid generation packages, such as Gridgen™, Grid Pro, and 
so forth, or use Laura’s interactive self .start utility to generate a 
single-block structured grid for simple families of 2D, axisymmetric 
and 3D blunt bodies — see Section 7.8 on page 59. 

Step la. Using self .start. To use self .start to generate a single- 
block grid, simply execute this interactive utility, e.g., 

°/ 0 self_start 

and answer all the questions. After a successful execution, 
this utility will have generated the following files: 

assign.tasks laura_bound_data 

laura.g laura_namelist_data self_start.log 

Examine the grid, laura.g, and proceed to Step 2. 

Step lb. External Grid Generation. Generate a single- or multi- 
block structured grid with the following rules: 

i. Right-handed grid coordinates 

ii. Longitudinal axis of the body aligned with the x-axis, 
oriented nose-to-tail 

and write the grid coordinates into a Plot3D file, laura.g. 
Run the interactive bounds utility (see Section 7.1 on page 57) 
and answer all the questions regarding the grid block topol- 
ogy: 

°/ 0 bounds 

This utility will automatically generate laura.bound.data, 
the connectivity file. 

Step 2. If you did not use self .start, create assign.tasks (see Section 7.7 
on page 59 for more info) and a laura.namelist.data file or copy 
the sample file from the [install.pref ix] /share/laura directory, 
where [install.pref ix] is the installation prefix specified when 
Laura was installed. Edit this file for your case — see Section 5.4 
on page 17 for more detail. 

8 See Section 8 on page 61 for complete worked examples and Appendix A on page 72 for 
how to restart cases run with versions of Laura prior to version 5. 
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Step 3. Create a tdata file (see Section 5.5 on page 35) to define the gas 
model condition for your specific simulation. 9 

Step 4. Run Laura, 

*/. laura 


or 


*/. mpirun -np [#] laura_mpi 

where # is the number of available processors. By the end of this 
step, the following hies will have been generated: 


laura_conv . out 
laura_new.g 
laura_surf ace . q 
laura_surf ace .nam 


laura.g.fvbnd 
laura_new.rst 
laura. q 
laura. nam 


Examine these hies before proceeding to the next step. 10 

Step 5. Change irest hag in the laura_namelist_data (see Section 5.4 on 
page 17) from 0 to 1, and copy the new generated grid and solution 
hies to laura.g and laura. rst hies; i.e., 

% cp laura_new.g laura.g 
*/. cp laura_new.rst laura. rst 

Step 6. Repeat the previous two steps until iterative convergence. 


9 A sample tdata is available in the [install_pref ix] /share/physics_modules instal- 
lation directory. The other datafiles that reside in this directory, e.g., kinetics_data, 
species_thermo_data, species_transp_data, and species_transp_data_0, may also be 
copied and tailored to suit a different thermodynamic model, curve-fit data, or thermo- 
chemical reactions are needed. See Section 5.5 on page 35 for more detail. 

10 See Section 6 on page 52 for complete description of laura output files. 


Users-Manual. tex 


11 


Rev: 43231, November 1, 2009 



5 Input Files 


Nominally, LAURA requires five input files as shown in the upper section of 
Table 1. Depending on the simulation requirements, however, other files may 
also be necessary and are shown in the second section of Table 1. All files are 
plain ASCII text unless otherwise noted. 

Table 1: 

: Laura input files. 

Filename 

Content 

assign.tasks 

laura.g* 

laura_bound_data 

laura_namelist_data 

tdata 

Sweep and relaxation directions 
PLOT3D grid 

Grid block face boundary conditions 
Simulation configuration 
Gas model 

hara_namelist_data 

kinetic_data 

laura.rstd 

laura.trn 

lauravtraj ectory_data 

laura_vis_data 

species_thermo_data 

species_transp_data 

species_transp_dataO 

surf ace_property_data 

Radiation mechanisms 
Specie reactants and products 
Flowficld solution for restart 
Transition location and length 
Trajectory points 
Viscous term treatment 
Specie thermodynamics 
Collision cross-sections 
High-order collision cross-sections 
Thermochemical surface properties 


Fortran unformatted binary, 3-D whole, multiblock Plot 3D. 
t Fortran unformatted binary. 


The following subsections describe all input files in detail, beginning with 
the nominally required files and then proceeding alphabetically as shown in 
the table. 

5.1 assign_tasks 

This file defines sweep and relaxation directions for each grid block. Each line, 
corresponding to each grid block, has six integers that are separated by at 
least one space. These integers correspond to nbk, mbk, mbr, lstrt, lstop, 
and mapcpu where 

nbk 

Block number. 
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mbk 

Sweep direction. The assigned value can be either 1, 2, or 3, correspond- 
ing to i-, j-, or /e-direction, respectively. 

mbr 

The line relaxation direction. Note: must be different than the sweep 
direction. 

Options are: 

0: Point-implicit, i.e., no line-relaxation 
1: Line-implicit along i coordinate 
2: Line-implicit along j coordinate 
3: Line-implicit along k coordinate 

lstrt 

The starting grid index in the sweep direction. Typically 1 . 

lstop 

The ending grid index in the sweep direction. Typically 0, which is 
shorthand for the maximum index. 

mapcpu 

CPU number for block nbk. Must be greater than zero. To run with a 
single processor, set this value to 1 for all grid blocks. 

When starting a new simulation where the /^-coordinate runs from the 
vehicle surface to the freestream boundary, sweeping in the /c-coordinate and 
solving point-implicitly is recommended, i.e. mbk = 3 and mbr = 0. After 
the shock has stabilized, switch to streamwise sweeps and solve line-implicitly 
along the ^-coordinate, i.e. mbk = 1 and mbr = 3. 

5.1.1 Example 1: Multiple Blocks per CPU or Vice-versa 

Suppose the grid has 2 blocks, but the number of available processors is not 
the same as the number of grid blocks. The user does not need to change any 
of the input files and can simply specify the number of processors available, 

e-g., 

I mpirun -np [#] laura_mpi 

where # is the number of available processors. Laura automatically assigns 
processors to blocks such that each processor receives approximately same 
number of grid cells. 
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5.1.2 Example 2: Deactivating Grid Blocks 


Suppose the grid has 50 blocks, but only blocks 20-25 need to be updated. In 
this case, assignxtasks would contain only the active blocks, e.g., 


20 3 0 1 

21 3 0 1 

22 3 0 1 

23 3 0 1 

24 3 0 1 

25 3 0 1 


0 1 

0 2 

0 3 

0 4 

0 5 

0 6 


nbk mbk mbr lstrt lstop mapcpu 


5.2 laura.g 

This hie is a multi-block, 3D-whole Plot 3D hie in Fortran unformatted binary 
format with double-precision reals. For convenience, here is a sample of the 
Fortran 95 code that Laura uses to read this hie: 

open ( 25, file=' laura.g 1 , f orm= ' unformatted 1 ) 
read(25) nblocks 

...allocate i,j,kblk(nb) and grid memory... 
read(25) (iblk(nb) , jblk(nb) ,kblk(nb) ,nb=l , nblocks) 
do nb = 1 , nblocks 

ixl = iblk(nb) ; jxl = jblk(nb) ; kxl = kblk(nb) 

. . .allocate grid(nb)%x,y,z memory. . . 

read(25) ( ( (grid(nb)%x(i , j ,k) , i=l , ixl) , j=l , jxl) ,k=l ,kxl) , & 

( ( (grid(nb)°/„y (i , j ,k) ,i=l,ixl) , j=l, jxl) ,k=l,kxl) , & 

( ( (grid(nb)%z(i, j ,k) ,i=l,ixl) , j =1 , jxl) ,k=l,kxl) 

end do 

A hie of this format, but named laura_new.g, n is generated by Laura at the 
end of a successful run. This hie is required and must have a right-handed 
coordinate system. Figure 1 shows the default laura coordinate orientation. 
Laura does not require a specihc coordinate or grid orientation but the angle- 
of-attack definition is predefined — see Section 5.4.3 on page 22 for more info. 


5.3 lauraJbouncLdata 

Grid block face boundary types are defined in laura_bound_data where each 
line corresponds to a grid block and contains six integers, one for each of the 
six faces: i min , i max , j min , j max , k min , and k max . Each integer specihes either 

n When running a trajectory sequence, the file will be named laura_####.g where #### 
is the trajectory point index. 
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Figure 1: Default LAURA coordinate system orientation. 


a physical boundary condition or block-to-block interfaces. An illustrative 
example is analyzed toward the end of this section. 

This hie is required and is generated automatically for grids created by 
Laura’s self _start utility — see Section 7.8 on page 59. This hie can also be 
created by using Laura’s interactive utility, bounds, by answering questions 
for each block. 

Valid face types are as follows: 


-9,...,0: Solid surface boundary. Up to ten different solid surface boundaries 
may be specihed. Thermochemical properties of solid surfaces that are 
different than type 0, which are specihed in laura_namelist_data, are 
defined in surf ace_property_data hie — see Section 5.15 on page 49. 

1: Outhow boundary (extrapolation). 

2: Symmetry boundary across y = constant. 

3: Farheld/Freestream boundary. 


4: Symmetry boundary across x = constant or z = constant. 

5: Reflection boundary across j — 1 symmetry (valid for axi-symmetric 
and/or 2D grids). 

6: Venting boundary. (See Section 5.4.16 on page 34 for more details.) 


7: Reflection boundary across i face singularity with periodic j boundary. 
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>1000000: This seven digit boundary number defines block-to-block face connectiv- 
ity. The first digit is always 1. The next three digits identifies the block 
number that is shared with the current block. The 5 th digit defines which 
i, j, or k face of the neighboring block is shared where 1 corresponds to 
imini 2 corresponds to i m ax, and so forth. The last two digits identify the 
relation of the remaining two indices: The 6 th digit can be either 1, 2, 
3, or 4 where values of 1 or 3 mean the first index of the host face is in 
the same direction as the first or second index of the neighboring face, 
respectively, and values of 2 or 4 mean the adjoining indices are in the 
opposite direction. The last digit can be either a 1 or 2 and indicates 
whether the second indices of the host and neighboring faces are in the 
same direction or they are in the opposite direction, respectively. 

For example, consider a block with the following laura_bound_data bound- 
ary condition numbers: 

1006421 1002111 2 1005121 0 3 

The first integer, 1006421, shows that i min of this block is shared with j max 
(fifth digit) of block 6 (the first three digits after 1). The first and second 
indices of the host face are j and k, respectively. Because the j index is 
connected to i, the first and second indices of the neighboring face are % and 

k, respectively, The 2 in the 6 th digit shows that the j index of the host face 
is in opposite direction of the i index of the neighboring face. The last digit, 

l, indicates that k indices of the host and and neighboring faces are along the 
same direction. This configuration is illustrated in Figure 2. 



Figure 2: Illustration used for block connectivity example. 

The second and fourth boundary-type integers can be explained similarly. 
The third boundary-type integer (corresponding to the j nan face) specifies a 
y-constant symmetry plane; the fifth boundary-type integer (corresponding to 
the k min face) specifies is a solid-wall boundary condition; and the last digit 
(corresponding to the k ma x face) specifies a freestream boundary. 
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5.4 laura_namelist_data 


Simulation configuration is specified through laura_namelist_data and is re- 
quired. This file is read as a Fortran 95 namelist and has the form, 

&laura_namelist 
velocity_ref = 5000.0 
density_ref = 0.023 
tref = 250.0 

alpha = 25.0 

[ variable = value 

/ 

where variable and their possible values are described in the following sec- 
tions, which are grouped according to farfield/freestream and aerodynamic 
coefficient reference quantities; thermochemical nonequilibrium flags; molec- 
ular transport flags; turbulent transport models flags; numerical parameters; 
grid adaptation, alignment, and doubling parameters; radiation and ablation 
flags; grid-file description; venting boundary condition flags; and solid sur- 
face boundary condition flags. Note that for all but the parameters shown 
in the above example, reasonable defaults have been chosen and only those 
parameters that differ from the defaults need to be specified. 

Detailed description of parameters and/or flags with their units and default 
values is presented under each of the above categories. The order of these 
parameters is arbitrary but is given here alphabetically for better readability. 

5.4.1 Ablation Flags 
abl at i on_opt i on 

An integer that specifies whether the pyrolysis ablation rate and wall 
temperature are computed in addition to the char ablation rate. This 
option only effects cases with bprime.f lag equal to 0 or 1. 

Options are: 

0: The pyrolysis ablation rate and wall temperature are computed, in 
addition to the char ablation rate, assuming steady-state ablation. 

1: The pyrolysis ablation rate and wall temperature are held constant 
(they are set to the values present in ablation_from_laura.m) while 
the char ablation rate is computed. 

ablation_verbose 

A logical flag to print out developer focused info on convergence of ab- 
lation. Default: .true. 


Free stream velocity, m/s 
reference density, kg/s 
Free stream temperature, K 
Angle-of-attack in xz plane, degrees 
Optional comment ] 


Users-Manual. tex 


17 


Rev: 43231, November 1, 2009 



blowing_model_0 

Character indicator for ablation model specification. Default: ‘ none ’ 
Options are: 

‘FIAT’ 

This ‘FIAT’ 12 option computes the blowing rate and surface tem- 
perature as a function of heating rate, pressure, and ablator ele- 
mental mass fractions. The user must specify the elemental mass 
fractions for char and pyrolysis gas or default of 100% carbon will 
be employed. 

‘ none 5 

No ablation-flowfield coupling. 

‘ specif ied’ 

This option specifies a blowing or suction rate as a function of pres- 
sure (see mdot_pressure_0). The user must specify the elemental 
mass fractions for pyrolysis gas or default of 100% carbon will be 
employed. Any specification for elemental mass fraction of char is 
ignored in this option. 

‘ equil_char_quasi_steady ’ 

This option solves the equilibrium surface ablation problem. The 
energy balance, elemental mass balance, and char equilibrium con- 
straint are solved to obtain the char ablation rate, wall temperature, 
and elemental composition at the surface. Along with the pressure 
from the normal momentum equation, these values define the equi- 
librium species composition at the wall. The user must specify the 
surface temperature type to be surface energy balance — see 
Section 5.4.10 on page 28 for more info. 

Presently, this model does not include an in-depth material response 
computation, which would provide the pyrolysis ablation rate and 
conductive heat flux into the solid material. These two values are 
required for solving the previously mentioned surface equations. To 
approximate these values, the steady-state ablation assumption is 
made, which specifies that the pyrolysis ablation rate is propor- 
tional to the char ablation rate and the in-depth conduction is pro- 
portional to the enthalpy at the surface. 

In the present framework of steady-state ablation, the ablation 
material is completely defined by CH0NSi_f rac_pyrolysis_0 (thus 

12 FIAT is a stand alone software and is needed if blowing_model_0 =< FIAT 1 . Request for 
access to FIAT can be made to NASA Ames Research Center. 
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CHONSi_f rac_char_0 is ignored in this option). These are defined 
below. Note that the computed ablation rate is the sum of the py- 
rolysis and char components. (See Sections 8.3 on page 67 and 8.4 
on page 69 for recommended procedure for an unspecified ablation 
computation.) 

‘ quasi_steady ’ 

This option specifies a quasi-steady ablation rate as a function of 
local pressure, heating rate, and temperature. The sublimation 
temperature and heat of ablation are specified by the user for a 
given ablator as a function of pressure. (See t_sublimation_0 and 
h_ablation_0.) If the surface temperature is less than the sub- 
limation temperature a zero blowing rate is defined. Otherwise 
the blowing rate is given by m = q/AH a M with appropriate non- 
dimensionalization employed before use. The user must specify the 
elemental mass fractions for pyrolysis gas or default of 100% car- 
bon will be employed. Any specification for elemental mass frac- 
tion of char is ignored in this option. The user must specify the 
surface temperature type to be surface energy balance — see 
Section 5.4.10 on page 28 for more info. 

bprime_f lag 

An integer defining if the b-prime approach is applied. Applicable only 
for blowing_model_0 = ‘ equil_char_quasi_steady J . See Section 5.4.1 
on the preceding page for more details of its application. Default: 1 

Options are: 

0: Do not use bprime approach, and instead use a rigorous diffusion 
model. This option is consistent with the “Fully-Coupled” approach 
defined in Ref. [2]. 

1: Use b-prime approach. This option is consistent with the “Partially- 
Coupled” approach defined in Ref. [2], 

2: Hold the ablation rate and wall temperature constant from the 
restart file, while applying the rigorous diffusion model (thus, the 
surface energy balance and char equilibrium constraint are not sat- 
isfied). This option is sometimes useful when transitioning from 
a bprirae.flag = 1 computation to a bprime_flag = 0 computa- 
tion. 

char_density_0 

Density of the char, kg/m 3 . Default: 256.29536, which is for the heritage 
AVCOAT. 
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CHONSi_f rac_char_0 


This rank 1 vector of extent 5 sets elemental mass fraction of C, H, 0, 
N, and Si species from char. Elemental mass fractions must be in this 
order and the sum of elemental mass fractions must equal 1.0 . Default: 
CH0NSi_f rac_char_0 = 1.0, 0.0, 0.0, 0.0, 0.0 

CH0NSi_f rac_pyrolysis_0 

This rank 1 vector of extent 5 sets the elemental mass fractions of pyrol- 
ysis gas species, which are C, H, O, N, and Si. Elemental mass fractions 
must be in this order and the sum of elemental mass fractions must be 1. 
Default is Graphite: 

CH0NSi_f rac_pyrolysis_0 = 1.0, 0.0, 0.0, 0.0, 0.0 
compute_mdot .initial 

An integer defining if the ablation rates are computed before the first 
flowfield iteration. 

Options are: 


0: Applies the ablation rates and wall temperatures present in the 

ablation_frora_laura.m hie. 

1: Computes the ablation rates and wall temperatures before the first 
flowfield iteration. 

freq.wall 

For bprime.flag = 1, it is an integer defining frequency of update to 
ablation wall boundary conditions. For bprime_flag = 0, it is an inte- 
ger defining frequency of update to the surface energy balance solution, 
which defines the wall temperature, while nexch defines the frequency of 
update to the remaining surface equations — see Section 5.4.8 on page 26 
for more info on nexch. Default: 50 
h_ablation_0 


A rank 1 vector of extent 3 used to compute the heat of ablation in 
MJ/kg for quasi_steady blowing option as 

h_ablation_0(l) + (h_ablation_0 (2) ) logp,,, , , 

+ (h_ablation_0 (3) )(logp u ,) 2 

where p w is the local pressure, in atmospheres. Example A H a u = 28.0 — 
1.375 log p w + 27. 2(logp UI ) 2 . Default: 0.0 


mdot_pressure_0 


A rank 1 vector of extent 2 is used to set the blowing or suction distri- 
bution defined as 


mdot_pressure_0(l) + (mdot_pressure_0(2)) 


P 


PooCJ 


( 2 ) 
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where p is the local pressure, p ^ is the reference density, and is the 
reference velocity. Positive value produces blowing distribution, while 
negative value produces suction distribution. Default: 0.0 

t_sublimation_0 

A rank 1 vector of extent 3 used to compute the sublimation temperature 
in degrees Kelvin for quasi .steady blowing option as 

t_subliraation_0(l) + (t_sublimation_0(2)) logjv, , . 

+ (t_sublimation_0(3))(logp u ,) 2 

where p w is the local pressure, in atmospheres. Example T su b = 3797.0 + 
342.0 log p w + 30. 0(logp M ,) 2 . Default: 0.0 

uncoupled_ablation_f lag 

An integer defining if an uncoupled ablation analysis is applied. The 
uncoupled ablation option is included to provide a baseline solution for 
the coupled ablation analysis. Default: 0 

Options are: 

0: Do not apply an uncoupled ablation analysis, ft means that the 
coupled ablation analysis discussed in Section 8.3 on page 67 is 
applied. 

1: Apply an uncoupled ablation analysis to a converged non-ablating 
flowheld. The procedure for applying the uncoupled ablation anal- 
ysis is discussed in Section 8.4 on page 69. 

virgin_density_0 

Density of virgin material, kg/m 3 . Default: 544.627742, which is for the 
heritage AVCOAT. 

5.4.2 Aerodynamic Coefficient Reference Quantities 
bref 

Yaw moment coefficient reference length, grid-units. Default: 1.0 
cref 

Pitching moment coefficient reference length, grid-units. Default: 1.0 
sref 

Reference area for aerodynamic coefficients, grid-units. Default: 1.0 
xmc 

^-coordinate of moment center, grid-units. Default: 0.0 
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ymc 


^/-coordinate of moment center, grid-units. Default: 0.0 
zmc 

z-coordinate of moment center, grid-units. Default: 0.0 

5.4.3 Farfield/Freestream Reference Quantities 
alpha 

Angle-of-attack in xz plane, degrees, such that 

u = cos (a) cos (yaw)] v = —sin(yaw)]w = sin(a)cos(yaw ) (4) 

where u, v, and w are velocities in x-, y-, and z-coordinate, respectively. 
Default: 0.0 

density_ref 

Free stream density, kg/m 3 . Default: 0.001 
rpm 

Spin rate, RPM. This is applicable only to axisymmetric cases. 

Default: 0.0 

tref 

Free stream temperature, K. Default: 200.0 
velocity_ref 

Free stream velocity, m/s. Default: 5000.0 
yaw 

Sideslip angle in xy plane, degrees. Default: 0.0 


5.4.4 Grid Adaptation, Alignment, and Doubling Parameters 
beta_grd 

This parameter controls grid points normal to the body (A;-grid points) 
are controlled by the following grid stretching function: 


K = i 


1 kmax 

3- 0+1 


- 1 


kmax 

P 1 kmax I 1 

P+1 1 


( 5 ) 


where k* are new grid points. The beta_grd parameter defines the (3 
coefficient in equation 5. No stretching will be performed if f3 < 1. 
Recommended value is 1.15, if used. Default: 0.0 
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epCLgrd 

Grid clustering around the shock is designed using the following function: 

k**i = e 0 k* 2 (l — kf)(k* + f sh) + k* (6) 

where k t refers to normalized value and k* is defined by equation 5 on the 
facing page. epCLgrd assigns the eo coefficient in equation 6. Maximum 
recommended value is 25/4, if used. Default: 0.0 

kmax_f inal 

The target number of grid cells in the h-direction. Triggered by the 
global L2 error norm set by kmax.error, cells along the k direction are 
increased by a factor of kmax.factor until reaching kmax_f inal. Any 
value less than the number of k grid cells in laura.g hie will be ignored, 
i.e., no coarsening. This option requires all blocks to be active. Default: 
0.0 

kmax_error 

When the global L2 error norm reaches this value, the number of grid 
cells in the k direction increases by a factor of kmax_f actor until the 
maximum allowable grid cells, kmax_f inal, is reached. This option re- 
quires all blocks to be active. Default: 0.01 

f ctr jmp 

This parameter is used to detect bow shocks. Bow shock is first detected 
when the sensing parameter defined by jumpflag exceeds by this value, 
while searching from inflow boundary. Default: 1.05 

frac_line_implicit 

This positive parameter, which must be less than or equal 1, sets a 
fraction of the line-implicit direction that is assigned in assign_tasks 
(see Section 5.1 on page 12). This parameter, which supersedes the given 
assignments in assignvtasks, will be applied to all the active blocks. 
This parameter is recommended where there might be an instability 
issue such as across strong shocks. Default: 0.7 

f sh 

Fraction of arc length distance between body and inflow boundary where 
bow shock is captured. Default: 0.8 

f str 

This parameter approximately defines fraction of Avcells in boundary 
layer region. Default: 0.75 
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jumpflag 


An integer flag to select the sensing parameter to be used in detecting 
position of captured shock for grid adaption. Default: 2 

Options are: 

0: Redistributes grid points in ^-direction for the target cell Reynolds 
number defined by re_cell parameter, without changing the do- 
main boundary. 

1: Selects pressure as the sensing parameter. 

2: Selects density as the sensing parameter. 

3: Selects temperature as the sensing parameter. 

4: Scales the grid distance in the ^-direction up by the value defined 
by fctrjmp parameter. 

maxmoves 

An integer number to assign maximum number of times that grid adap- 
tion is performed. The value of zero, however, can be used for unlimited 
number of times. Default: 0 

max_distance 

This real number defines maximum distance from the body surface that 
grid outer boundary can be moved away by any one of the flags. This is 
often useful especially when adapting shock into wake, where the adapt- 
ing grid may become skewed due to presence of sharp gradients. This 
value defines the maximum length of the wake. Default: 1.0E+6 

movegrd 

An integer number for number of cycles between each grid alignment. A 
zero value disables any grid alignment. Default: 0 

re_cell 

This value defines the target cell Reynolds number at the wall after a 
grid movement. Note : If the grid is moving radically Default: 0.1 Note: 
Use a higher value for re.cell for the first grid alignment, if the grid 
movement causes radical changes in the grid. 

The cell Reynolds number is defined as 


Recell = pcAn/ n (7) 

Here p is the flow density, c is the sound speed, An is the cell height, 
and // is the flow viscosity. 
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5.4.5 Grid File Description 
dimensionality 

A string flag to select the dimensions of the problem. Default: ‘3DL 
Available options are: 

‘ axisymmetric J 

This option selects axisymmetric flow solution. This requires a 
domain with single-cell width in the j-direction. 

‘ 2D ’ 

This selects two-dimensional problem, which requires a domain with 
single-cell width in the j-direction. 

‘ 3D ’ 

This option solves three-dimensional problems. 
grid_conversion_f actor 

This parameter scales the grid to meter. One grid unit equals grid_conversion_f actor 
meters. Default: 1 

single_precision_grid 

Logical flag: .true, if grid points in the laura.g hie are stored as single 
precision values, otherwise double precision is assumed. Default: .false. 

5.4.6 Initialization 
init_vel_f ctr 

A real number between 0 and 1 to reduce the initial velocities u, v, w 
within the domain to avoid creating a vacuum for wake how problems, 
which may lead to an invalid solutions. The initial near zero velocity has 
also shown to ease up the formation of bow shock. Default: 0.01 

5.4.7 Molecular Transport Flags 
ivisc 

An integer hag to engage viscous terms (inviscid=0/viscous=2). Default: 

2 

mass_driven_dif f 

A logical hag to engage binary diffusion driven by mass fraction gradient. 

Default: .false. 

multi_component_dif f 

A logical hag to engage multi-component diffusion by Stefan-Maxwell 
equation sub-iterations. Default: .false. 
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mole_driven_dif f 

A logical flag to engage binary diffusion driven by mole fraction gradient. 
Default: .true. 

navier_stokes 

A logical flag to select equation set for Thin-Layer Navier Stokes or Full 
Navier Stokes. The navier_stokes = .false, may be used to select 
Thin-Layer Navier Stokes. Default: .false. 

pressure.dif fusion 

A logical flag to engage pressure diffusion term on Stefan-Maxwell ap- 
proximation. Default: .false. 

schmidt _number 

A constant Schmidt number may be specified to calculate diffusivities. 

If the value is negative, diffusivities are computed directly from collision 
cross sections. Default: -1.0 

prandtl_number 

A constant Prandtl number may be specified to calculate conductivi- 
ties. If the value is negative, conductivities are computed directly from 
collision cross sections. Default: -1.0 

5.4.8 Numerical Parameters 
cfll 

Initial value of CFL number. Default: 5.0 
cf 12 

Final value of CFL number. Default: 5.0 
epsa 

Eigenvalue limiting factor. Default: 0.3 
hrs 

The maximum total CPLT time for simulation, hours. Default: 10 
iramp 

Number of cycles to ramp from cfll to cfl2. Default: 200 
irest 

An integer flag to start the simulation either using freestream values 
(irest = 0), or using the existing solution from laura.rst file (irest = 1). 
Default 0 
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iterwrt 

Number of cycles between saves of all output files — see Section 6 on 
page 52. Default: 200 

j update 

Number of cycles between jacobian updates. Default: 10 
ncyc 

Number of global iterations. Default: 1000 
nexch 

Number of cycles between exchange of data between processors and up- 
dating boundary conditions. Default: 2 

nitf o 

Number of cycles using lst-order spatial accuracy. Default: 0 
nordbc 

Boundary condition calculation using l st -order (nordbc = 1) or 2 nd - 
order (nordbc = 2) accuracy. Default: 1 

ntran 

Number of cycles between transport properties update. Default 1 
rf _inv 

Inviscid relaxation factor. Default: 3.0 
rf _vis 

Viscous relaxation factor. Default: 1.0 
rf _chem 

Chemical source term reduction factor sometimes useful to ’’ease in” 
simulations very close to equilibrium. This factor, which must be greater 
than 1.0 when it is used, changes the answers and must ultimately equal 
1.0 in the final simulation. Default: 1.0 

rmstol 

A real number to stop the iterations after L2 norm residual reaches this 
value. Default: 1.0E-10 

sub_iterations 

An integer number to control multiple passes through the linear solver, 
which has been found to increase robustness for some high energy appli- 
cations. Default: 1 


Users-Manual. tex 


27 


Rev: 43231, November 1, 2009 



5.4.9 Radiation Flags 
radiation 

A logical flag to enable coupling between radiation equation(s) and flow 
equations. See Section 8.2 on page 66 for coupled radiation procedure. 
Default: .false. 

radiation_input_only 

A logical flag to enable creation of input file hara_out.m to compute 
radiation outside of LAURA when radiation is .false. Default: .false. 

maxrad 

An integer number to assign maximum number of calls to radiation in- 
terface. The value of zero can be used for unlimited number of times. 
Default: 0 

nrad 

Number of iterations between calls to Hara. Default: 3000 
iincurad, jincurad 

Increment between i and j lines, respectively at which radiation profile 
is computed. Interim lines are interpolated. Default: 3 

frac_rad_new 

Relaxation factor on radiation. ' V(q ra d ) = fracjrad-newV (q ra d) n + (1 — 
fracj-adjT,ew)'V(q ra d) n ~ 1 - Default: 1.0 

absorptivity 

Fraction of incident radiative energy absorbed by the wall. Affects sur- 
face energy balance and computation of radiative equilibrium wall tem- 
perature. Default: 1.0 

tw_rad_f lag 

An integer flag to engage Tauber- Wakefield approximation for radiation 
cooling on surface-energy balance (on=l/off=0). Default: 0 

5.4.10 Solid Surface Boundary Condition Flags 
catalysis unodelJ) 

A character identifier for selecting catalysis model for surface type 0 (see 
Section 5.3 on page 14). This flag is good only for multi-species reacting 
gases and will be ignored for single-specie gas models. The catalysis 
model name must be surrounded by quotation marks, e.g., ‘ ’. Default: 
‘ super-catalytic ’ 

Available catalysis models are: 
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' CCAT-ACC ’ 

This option uses the following surface catalytic coefficients [3] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 


f 13.5e- 835 °/ T - T* < 1359.0 K 

\ 5.0 x 10“ 8 e 18023 / T ™ T* > 1359.0 K 


_ f 4.0e -7625 / T ™ T* < 1475.0 K 

lN ~ \ 6.2 x I0“ 6 e 1210 °/ T » T* > 1475.0 K 

where T* is calculated as 


f min( max(1255.0, T w ), 1659.0) for 
| min( max(1255.0 ,T W ), 1900.0) for 7 ^ 


(9) 


( 10 ) 


‘ CSiC’ 

This option uses the following surface catalytic coefficients [3] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 

70 = 6.415 x io- 4 e 3498 - 4/T » (11) 

= 3.993 x 10” 4 e 4402/ ' 7 ^ (12) 

where T* is given as 

T* — min( max(1100.0, T w ), 1920.0) (13) 


' CSiC-SNECMA ’ 

This option uses the following surface catalytic coefficients [3] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 


1o = 1n = 9.593 x i(r 5 e 7002 - 9/T “ (14) 

where T* is given as 

T* = min( max(1350.0, T w ), 1920.0) (15) 

‘ equilibrium-catalytic ' 

This option sets species concentrations to equilibrium values at wall 
pressure and temperature based on elemental mass fractions in the 
cell above the solid surface boundary. 
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‘ f ully-catalytic ’ 

This option assumes that all the atomic and ionized oxygen, nitro- 
gen, carbon, and so forth catalyzes to molecular oxygen, nitrogen, 
carbon, and so on, respectively; i.e., 


7o = 7n = 7c = • • • = 1 


(16) 


‘ non-catalytic ’ 

This option assumes that no atomic or ionized oxygen, nitrogen, 
carbon, and so forth catalyzes to molecular oxygen, nitrogen, car- 
bon, and so on, respectively; i.e., 


7o = 7n = 7c = ••• = 0 


( 17 ) 


‘ RCC-LVP ’ 

This option uses the following surface catalytic coefficients [3] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 


7 o 


7 N 


7.5e -8283 / T ™ 

T* < 1499.0 K 

(18) 

2.5 x io-V 7533 / T » 

T* > 1499.0 K 

6.0 x 10~ 2 e^ 2605 / T ^ 

T* < 1529.0 K 

(19) 

1.5 x 10 -5 e 1008 °/ T ™ 

T* > 1529.0 K 


where T* is calculated as 


rjn* 


min( max(1255.0, T w ), 1799.0) for 70 
min( max(1255.0, T w ), 1954.0) for 7^ 


( 20 ) 


' Scott-RCG’ 

This option uses the following surface catalytic coefficients [4] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 

70 = 16. Oe -10271 / 7 ™ 1400 <T W < 1650 (21) 

7^ = 7.14 x i 0 - 2 e- 2219 -°/^ 1090 < T w < 1670 (22) 

The same equations will be used even if the wall temperature, T w , 
is out of the specified range, in which case a warning will be issued 
to the stdout. 

‘ Stewart-RCG’ 

This option uses the following surface catalytic coefficients [3] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
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nitrogen, respectively. 


7 o 


7 N 


5.0 x 10 3 e 40 °/ T ™ 
1.6 x io~ 4 e 1326 /' ru ' 
5 2e -8835 / T “ 

39 x iO” 9 e 2141 °/ T ™ 

5.0 x 10 -4 

2.0 x lO“ 5 e 150 °/ T “ 
10 Oe -10360 / T “ 

6.2 x 10 _6 e 1210 °/ T ™ 


T w < 502 
502 < T w < 978 
978 < T w < 1617 
1617 < T w 

T w < 465 
465 <T W < 905 
905 <T W < 1575 
1575 < T w 


(23) 


(24) 


‘ super-catalytic 5 

This option sets the species mass fractions to free stream values as 
defined in tdata — see Section 5.5 on page 35. 

‘ Zoby-RCG ’ 

This option uses the following surface catalytic coefficients [5] for 
catalyzing atomic oxygen and nitrogen to molecular oxygen and 
nitrogen, respectively. 

7 o = 9.41 x io- 3 e - 658 - 9 / T » 900 <T W < 1500 (25) 

7iv = 7.14 x i 0 - 2 e- 2219 -°/^ 1090 <T W < 1670 (26) 

The same equations will be used even if the wall temperature, T u , , 
is out of the specified range, in which case a warning will be issued 
to the screen. 

emiss_a_0, emiss_b_0, emiss_c_0, emiss_d_0 

Real number values to calculate cmissivity, e, for solid surface boundary 
type (see Section 5.3 on page 14) from the following equation: 

e = e a + tbT w + e cT^ + qT 3 (27) 

where T w is the surface temperature. Values for e a , e&, e c , and are de- 
fined by emiss_a_0, emiss_b_0, emiss_c_0, and emiss_d_0, respectively. 
Default: emiss_a_0=0.89, emiss_b_0=0.0, emiss_c_0=0.0, emiss_d_0=0.0 

ept 

Under- relaxation parameter for radiative equilibrium wall temperature: 

CP = (! “ e P + e P tQij 1 ( 28 ) 

where n denotes iteration level, g" 4 " 1 is the most recent value of convective 
heat flux, q, h3 , qfj is the value of convective heat flux from the previous 
iteration as adjusted by previous application of this formula, and qfj 1 is 
the newly adjusted value for convective heat flux at the n+ 1 iteration 
level. Default: 0.01 
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surf ace_temperature_type_0 


Character identifier for surface temperature model selection. Default: 

‘ constant ’ 

Options are: 

‘ adiabatic ’ 

Surface temperature will be such that there is no heat transfer be- 
tween the surface and the gas adjacent to the surface. 

‘ constant ’ 

The surface temperature stays constant as given by twallibc value, 
‘radiative equilibrium ’ 

The surface temperature is calculated so that the heat flux to the 
wall, q w , is in equilibrium with radiation heat flux: 

q w = eaT, * (29) 

where a is the Stefan-Boltzmann constant, and e is the surface 
emissivity defined by emiss_a_0, emiss_b_0, emiss_c_0, emiss_d_0 
values. 

‘surface energy balance’ 

This option is required for ablating surfaces to compute surface 
temperature as a function of the surface energy balance and the 
relevant surface chemistry kinetics. 

surf ace_group_name_0 

Character descriptor for surfaces with solid surface boundary types (see 
Section 5.3 on page 14). Any character can be specified to group solid 
surface boundaries. Default: ‘default surface 0’ 

twallibc 

Initial wall temperature for solid surface boundaries, K. (See Section 5.3 
on page 14.) The wall temperature stays constant as specified by this pa- 
rameter if surf ace_temperature_type_0 = ‘constant’. Default: 500.0 

5.4.11 Surface Recession Flags 
shape_change 

A logical flag engaging the geometry shape change due to ablation. De- 
fault : .false. A trajectory file is required if shape_change=.true. - 
see Section 5.10 on page 44. 
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5.4.12 Thermochemical Nonequilibrium Flags 
chem_f lag 

An integer flag to engage chemical source term for non-equilibrium flow 
(on=l/off=0). Default: 0 

therm_f lag 

An integer flag to engage thermal source term for non-equilibrium flow 
(on=l/off=0). Default: 0 

5.4.13 Time Accurate Flags 
itime 

An integer flag to engage (itime = 1) or disable (itime = 0) time ac- 
curate simulation. Default: 0 

subiters 

An integer number to specify number of iterations between each time 
step for time-accurate simulations. Default: 10 

5.4.14 Trajectory Related Flags 
traj ectory_data_point 

Pickup simulation from this line in the file laura_trajectory_data. 
Default: 0 if laura_trajectory_data is not present. 

Default: 1 if laura_trajectory_data is present. 

Note that the reference quantities are defined consistently across the 
trajectory using the values supplied in namelist above. The reference 
quantities are NOT reset according to the time varying free stream con- 
ditions. Consequently, dimensionless values of density and velocity in 
free stream will not generally equal 1. Dimensionless values of total en- 
thalpy in free stream will not equal 0.5 which may impact post-processing 
tools that are key on a specific number for total enthalpy to detect the 
boundary-layer edge. 

The algorithm is most efficient if the restart solution is converged (or 
nearly converged so that line-implicit relaxation may be applied) at free 
stream conditions equal to the reference quantities in laura_namelist_data. 
If one wishes to pickup a computation for trajectory_data_point > 1 
then it is best to start at the converged restart file for the previous 
trajectory point. Restart files and post processing files for each tra- 
jectory point have the trajectory point number included as part of the 
root name. Thus if one wants to pickup at trajectory point 12 one 
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should copy laura_0011.rst to laura.rst and (if grid is being updated) 
laura_0011.g to laura.g. This advice is offered because restart solutions 
are converted according to the ratio of density and velocity between ad- 
jacent trajectory points to bring interior initial conditions closer to the 
new inflow boundary conditions. 

5.4.15 Turbulent Transport Models 
prandtl_turb 

Turbulent Prandtl number. Default: 0.9 
schmidt_turb 

Turbulent Schmidt number. Default: 0.9 
turb_int_inf 

Farhcld turbulent intensity non-dimensionalized by square of the refer- 
ence velocity. Default: 0.001 

turb_model_type 

An integer flag to select turbulence modeling. Default: 0 
Options are: 

-1: Baldwin-Lomax algebraic turbulence model. [6] 

-2: Cebeci-Smith algebraic turbulence model. [6] 

0: Laminar flow, i.e. , no turbulence model. 
turb_vis_ratio_inf 

Farheld ratio of turbulent to laminar viscosity. Default: 0.01 
xtr 

Transition location along x-axis in grid units. Default: 0.0 

5.4.16 Venting Boundary Condition Flags 
vacuum_pressure_coef f icient 

This parameter, which is in N/m 2 sets the pressure coefficient behind 
type 6 boundaries (see Section 5.3 on page 14), which forces flow out of 
the domain. This coefficient is defined as 

vacuum_pressure_coeff icient = (30) 

where p 0 is back pressure where the gas is venting out, and p ^ is farheld 
pressure. Default: 0.0 
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vacuum_pressure_f actor 

Factor on pressure across type 6 boundary to force flow out of domain. 

vacuum_pressure_f actor = — , (31) 

Pi 

where p\ is the pressure just before the gas vents out. Default: 0.01 


5.5 tdata 

The gas model is defined in this hie. It contains a list of key words, sometimes 
followed by numeric values, which identify components of the gas model. One 
or more spaces must be present between keyword and values when appearing 
on the same line. Spaces may appear to the left or right of any key word. The 
first line of the hie must not be blank, however. 

The following subsections describe available gas model options. 


5.5.1 Perfect Gas 


The perfect-gas option is engaged with either of the following keywords: 
perfect_gas, PERFECT_GAS, Perfect.Gas, or Perfect_gas. 

If no further data is provided in this hie, this single line tdata hie will 
assume the following parameter values in SI units: 


gamma 

mol_wt 

sutherl 

suther2 

prand 


1.4 

28.8 

0. 1458205E-05 
110.333333 
0.72 


Here, gamma is the gas specihc heat ratio, mol_wt is the gas molecular weight, 
prand is the gas Prandtl number, and sutherl and suther2 are the hrst and 
second Sutherland’s viscosity coefficients, Si and s 2 , respectively, defined as 

T 3/2 


fl 5l 


T + s 2 


(32) 


These values can be modified and explicitly defined in the tdata by the 
keyword &species_properties in the second line followed by the gas param- 
eters and / at the last line of the hie. For example, 


perf ect_gas 
fespec ies_properti.es 
gamma = 1.4 
mol_wt = 28.0 
sutherl = 0.1E-05 
suther2 = 110.3 
prand = 0.7 
/ 
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5.5.2 Equilibrium Gas 

To engage the Tannehill curve fits for thermodynamic and transport properties 
of equilibrium air [7], the following keyword should be used in the first line of 
the tdata file: equilibrium_air_t. No additional inputs or files are required 
to engage the Tannehill option for equilibrium air. 

To use a table look-up capability for equilibrium gases [8], the following 
keyword should be placed in the tdata file, instead: equilibrium_air_r. Note 
that this option still uses the Tannehill transport properties. 

5.5.3 Mixture of Thermally Perfect Gases 

If the gas is a mixture of thermally perfect gases and multi-species transport 
solution is desired the species names followed by their mass fractions must be 
provided in the tdata file. Thermal state of the gas may be defined as the 
first entry by either of the following flags: 

one 

This flag assumes that all the species are thermally in equilibrium state. 
That is translational temperature, T, and vibrational temperature, T v 
are equal. This is known as one-temperature, 1-T, model. 

two 

This flag assumes that energy distribution in the translational and rota- 
tional modes of heavy particles (not electrons) are equilibrated at trans- 
lational temperature, T, and all other energy modes (vibrational, elec- 
tronic, electron translational) are equilibrated at vibrational tempera- 
ture, T v . This is known as two-temperature, 2-T, model. 

FEM 

This option, called Free-Energy Minimization, causes the species conti- 
nuity equations to be replaced with elemental continuity equations and 
equilibrium relations for remaining species. 

One temperature model is assumed if thermal state of the gas is not provided 
in the first line of the tdata file. In this case, the first line must contain species 
information. Note, the first line must not be blank. 

Subsequent file entries include species names and their mass fractions at 
freestream/farfield boundary. Only one specie per line is allowed. The species 
mass fraction at the boundary is defined in the same line as the species name 
separated by one or more spaces. If no value appears to the right of the species 
name then that species is assumed not to be present at the boundary but may 
be produced through chemical reactions elsewhere in the flowfield. 
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Example 1: 1-T, 5-species air model: In this example, only molecular 

oxygen and nitrogen are present on freestream/farfield boundary, but atomic 
nitrogen and oxygen and nitric oxide may be produced elsewhere in the flow 
field due to chemical reactions. 


one 


N2 

.767 

N 


02 

.233 

0 


NO 



Example 2: 2-T, 11-species air model: In this example, the gas is as- 

sumed to be a mixture of 11 thermally perfect gases. A solution to a thermal 
non-equilibrium state of the gas is also desired (2-T model). 

two 

N2 .767 
N 

02 .233 

0 

NO 

02 + 

0 + 

N0+ 

e- 

5.6 hara_namelist_data 

This file controls the radiation models used by the Hara radiation mod- 
ule [9,10]. It is optional for coupled radiation simulations. If it is not present, 
then the code automatically chooses the radiative mechanisms associated with 
species present in the flowfield (and have number densities greater than 1000 
particles/cm 2 ), and other options are set to the defaults listed in the following 
section. For users not experienced in shock-layer radiation, it is recommended 
that this hara_name list.dat a file not be applied (meaning it is removed from 
the working directory), therefore allowing the radiative mechanisms to be au- 
tomatically chosen and the default model options applied. 

5.6.1 Specifying radiation mechanisms for atomic species 

The treatment of radiation resulting from atomic lines, atomic bound-free 
and free- free photoionization (referred to here as atomic continuum), and neg- 
ative ion continuum is available for atomic carbon, hydrogen, oxygen, and 
nitrogen. These mechanisms are specified through the following binary flags 
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(on=l/off=0). If any of these flags are not present in hara_namelist_data, 
then that flag is set to true only if the number density of the associated atomic 
specie is greater than 1000 particles/cm 2 somewhere in the flowheld. 

treat, [?] -lines 

A binary flag to enable the treatment of atomic lines for specie [?] , 
where [?] can be c, h, n, and o, for atomic carbon, hydrogen, nitrogen 
and oxygen, respectively. 

treat, [?] ,cont 

A binary flag to enable the treatment of atomic bound-free and free-free 
continuum for specie [?] , where [?] can be c, h, n, and o, for atomic 
carbon, hydrogen, nitrogen and oxygen, respectively. 

treat, [?] _other 

A binary flag to enable the treatment of the negative-ion continuum for 
specie [?] , where [?] can be c, h, n, and o, for atomic carbon, hydrogen, 
nitrogen and oxygen, respectively. 

5.6.2 Specifying radiation mechanisms for molecular species 

The treatment of radiation resulting from numerous molecular band systems 
is available through the following binary flags (on=l/off=0). If any of these 
flags are not present in hara_namelist_data, then that flag is set to true 
only if the number density of the associated molecular specie is greater than 
1000 particles/cm 2 somewhere in the flowheld. Additional band systems are 
listed in Appendix B on page 75. These additional band systems are generally 
considered negligible relative to those listed in this section, and therefore for 
computational efficiency, they are not engaged by default. Definitions of each 
band system and the modeling data applied are discussed in Refs. [9,11]. 

treat _band_c2_swan 

A binary hag activating the C2 Swan band system. 

treat _band_c2h 

A binary hag activating the C2H band system. 

treat_band_c3 

A binary hag activating the C3 and Vacuum Ultra-Violet (VUV) band 
systems. 

treat_band_cn_red 

A binary hag activating the CN red band system. 
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treat _band_cn_violet 


A binary flag activating the 
treat _band_co4p 
A binary flag activating the 
treat_band_co_bx 
A binary flag activating the 
treat_band_co_cx 
A binary flag activating the 
treat _band_co_ex 
A binary flag activating the 
treat _bandJi2_lyman 
A binary flag activating the 
treat _band_h2_werner 
A binary flag activating the 
treat _band_n2fp 
A binary flag activating the 
treat _band_n2sp 
A binary flag activating the 
treat _band_n2pfn 
A binary flag activating the 
treat_band_n2_bhl 
A binary flag activating the 
treat_band_n2_bh2 
A binary flag activating the 
treat_band_no_beta 
A binary flag activating the 
treat_band_no_delta 
A binary flag activating the 
treat _band_no_epsilon 
A binary flag activating the 


CN violet band system. 
CO 4+ band system. 


CO B-X band system. 


CO C-X band system. 
CO E-X band system. 


H 2 Lyman band system. 


H 2 Werner band system. 
N 2 1+ band system. 


N 2 2+ band system. 

N^~ first- negative band system. 


N 2 Birge-Hopheld I band system. 


N 2 Birge-Hopheld II band system. 
NO beta band system. 


NO delta band system. 


NO epsilon band system. 
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5.6.3 Atomic line models 


There are various models available for atomic line radiation, one of which 
must be chosen for each specie that engages atomic line radiation (as specified 
using treat. [?] _lines). This choice of atomic line model is made using the 
following flags. The listed defaults are applied if the individual flag is not 
present in hara_namelist_data, or if hara_namelist_data is not present in 
the working directory. All model types in this category must be surrounded 
by a quotation marks, e.g. ‘ ’. 

c_atomic_line_model, h.atomicLLine .model 

A character identifier for selecting the atomic line model for atomic car- 
bon or hydrogen. Presently, the only available option is the model com- 
piled in Ref. [11], which is referred to here as the Complete Line Model 
(CLM). Default : 'elm' 

n_atomic_line_model, o_atomic_lineunodel 

A character identifier for selecting the atomic line model for atomic ni- 
trogen or oxygen. The available models are compiled and compared in 
Ref. [9], which is referred to here as the Complete Line Model (CLM). 
Default : ‘elm’ Available models are: 

‘all multiplets’ 

This model treats all lines as grouped multiplets. This significantly 
reduces the number of lines treated as well as the computational 
expense. However, this grouped multiplet approximation will lead 
to errors for non-optically-thin conditions. 

‘ elm’ 

This model, which stands for Complete Line Model, applies the 
individual treatment of strong atomic lines while applying multiplet 
averages for weak lines. This is the recommended model. 

5.6.4 Electronic state population models 

These flags specify the model applied for predicting the electronic state popu- 
lations of atoms and molecules. The listed defaults are applied if the individual 
flag is not present in hara_namelist.dat a, or if hara_namelist_data is not 
present in the working directory. All model types in this category must be 
surrounded by a quotation marks, e.g. ‘ ’. 

Atomic electronic states 

The electronic state populations for atoms are required for computing atomic 
line and photoionization emission and absorption. The compilation and com- 
parison of the available models are presented in Ref. [10]. 
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c_electronic_state, h_electronic_state 


A character identifier for selecting the electronic state model for atomic 
carbon and hydrogen. Available models are (Default : ‘boltzmann’): 

‘ boltzmann’ 

Applies Boltzmann population of electronic states. 

‘ Gally_lst_order_LTNE’ 

Applies the Gaily first-order local thermodynamic nonequilibrium 
method [13], which approximately accounts for the non-Boltzmann 
population of atomic states. 

n_electronic_state, o_electronic_state 

A character identifier for selecting the electronic state model for atomic 
nitrogen and oxygen. Available models are (Default : ‘CR’): 

‘ boltzmann’ 

Same as for c_electronic_state 

‘ Gally_lst_order_LTNE’ 

Same as for c_electronic_state 

‘ CR’ 

Applies the detailed Collisional Radiative (CR) model developed in 
Ref. [10], 

‘ AARC ’ 

Applies the Approximate Atomic Collisional Radiative (AARC) 
model developed in Ref. [10]. This model is essentially a curve-fit 
based approximation of the CR model, which allows for improved 
computational efficiency with a slight loss in accuracy. 

Molecular electronic states 

The electronic state populations for molecules are required for computing 
molecular band emission and absorption. The compilation and comparison 
of the available models are presented in Refs. [10,14], 

molecular_electronic_state 

A character identifier for selecting molecular electronic state for all molec- 
ular band systems. Available models are (Default : ‘CR’): 

‘ boltzmann’ 

Applies Boltzmann population of electronic states. 
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‘ CR 


Applies a detailed Collisional Radiative model considering both 
heavy-particle and electron impact transitions. Some molecular 
states are still assumed Boltzmann with this model because no data 
is presently available for the CR model. 

5.6.5 Other flags 
usevtriangles 

A logical flag specifying whether optically-thin atomic lines are modeled 
as triangles to reduce computational time. This option has shown to 
result in a negligible loss of accuracy while greatly reducing the compu- 
tational time, [9] and is therefore recommended. Default : .true. Note: 
This flag is automatically set to .true, when n_ or o_atomic_line_model= 
‘ elm’ — see Section 5.6.3 on page 40. 

use_edge_shift 


bound-free radiation. (on=l/off=0). Default : .true. 

5.7 kinetic_data 

This hie defines possible chemical reactions and is optional. By default, 
kinetic.data is read from [install-prefix] /share/physics_modules, but 
if present in the local run directory, LAURA will read reactions from the local 
hie instead . 13 

Reactants and products can be any species defined in the species_thermo_data 
hie — see Section 5.12 on page 46. A sample entry looks like this, 

02 + M <=> 20 + M i 

2 . 000e+21 -1.50 5.936e+04 2 

teffl =2 3 

expl =0.7 4 

t_eff_min = 1000. 5 

t_eff_max = 50000. 6 

C = 5.0 7 

0 = 5.0 s 

N = 5.0 9 

H = 5.0 10 

Si = 5.0 11 

e- = 0 . 12 

13 The precise installation location is given by Laura during startup. It can also 
be found on Unix-like systems from the executable itself by issuing the command 

strings laura I grep share. 
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The first line specifies the reaction while line 2 provides three coefficients of 
an Arrhenius- like equation, 


K, = c f T; fI e-- /kT ’» (33) 

where Cf is the pre-exponential factor, rj is the power of temperature depen- 
dence on the pre-exponential factor, eo is the Arrhenius activation energy, and 
k is the Boltzmann constant. The arrowheads in line 1 signify the allowed 
directionality of the reaction. The symbol => denotes forward reaction only 
while <=> denotes forward and backward rates are computed. The coeffi- 
cients in line 2 correspond to Cf, rj, and eo/h, respectively. For reactions with 
a generic collision partner, M, such as this one, these coefficients correspond 
to Argon; and other collision partners and their efficiencies (multipliers of c/) 
are specified on lines following line 5 and 6, which give the valid temperature 
range for the reaction. The effective temperature, T e // , is defined according 
to a given integer number in line 3; Default: 2 


teff l,tef f 2 

Flag defining formula to compute the effective temperature T e ff for the 
forward rate and backward rate, respectively. It is engaged for the case 
of thermal nonequilibrium. Options for teff are: 

1: Teff = T tr 

Q. 7^ rpexplrpl— expl 

L • 1 eff -Ttr 1 V 

3: Teff = T v 

where T tr and T v are translational and vibrational temperatures, respec- 
tively. Default: 1 

expl 

The exponent used to define the effective temperature when teffl = 2 
(forward rate) or teff 2 = 2 (backward rate). See previous equations for 
teff options. Default: 0.7 

t_ef f _min 

The minimum temperature for T e ff to compute reaction rates to circum- 
vent stiff source terms. Default: 1000. 

t_ef f _max 

The maximum temperature for T e ff to compute reaction rates to cir- 
cumvent stiff source terms. Default: 50000. 
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5.8 laura.rst 

This Fortran-unformatted binary file has the flowfield solution for every grid 
cell, boundary surface values, and grid cell derivatives for each block face. 
The number of variables in this hie varies depending on number of conserva- 
tion equations that LAURA needs to solve as specified through the tdata and 
laura_namelist_data hies. A hie of this format, but named laura_new.rst, 14 
is generated by LAURA at the end of a successful run. The laura.rst hie is re- 
quired only if restarting from a prior run, i.e., irest = 1 in laura_name list.dat a 
as described in Section 5.4.8 on page 26. 

5.9 laura.trn 

Turbulent transition location and transition length are specihed in laura.trn. 15 
One line of data consisting of four integers followed by at least one space with 
the following entries per block is required: block_number, turbulent_f lag, 
transition_location, and transition_length_f actor where 
turbulent f 1 ag 

One of -1, 0, or 1, which correspond to laminar, transition location in 
block, or fully turbulent how, respectively. 

transitiondLocation 

Specihed as an x coordinate. This value used only if turbulent.f lag 
is 0. 

trans it ion_length_f actor 

Defined as (transition length)/!/, where L is the distance from the nose 
to the transition location. Use 0 for instantaneous transition. 

5.10 laura_trajectory_data 

This hie is used to sequentially simulate points along a trajectory. Each line of 
the hie defines a single trajectory point. The trajectory point data is entered 
in free format with time, velocity, density, temperature, alpha, and yaw in 
MKS units (6 entries per line). The simulation will start at trajectory point 1 
by default unless another line number is specihed in the laura_namelist_data 
using the namelist variable trajectory_data_point — see Section 5.4.14 on 
page 33. 

The restart solution hies and post processing hies for each trajectory point 
are saved with a four digit number added to the usual root name of these hies 
corresponding to the trajectory point number. 

14 When running a trajectory sequence, the file will be named laura_####.g where #### 
is the trajectory point index. 

1o Laura transition files from versions prior to 5 can also be used. 
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5.11 laura_vis_data 


This file contains the data set that overrides the default viscous terms in the 
i-, j-, and /^-directions. There are four integers consisting of block number 
and toggles for each viscous term per line in this hie: blk_num, ivis, jvis, 
and kvis where ivis, jvis, and kvis are viscous terms in the i-, j-, and 
^-directions, respectively. The viscous-term flag options are: 

0: Viscous terms are not engaged in the respective direction; i.e. inviscid 
flow. 

1: Viscous terms and a reduced eigenvalue limiter are engaged in the respec- 
tive direction. The reduced eigenvalue limiter is to prevent distortion of 
computed heating. 

2: Viscous terms are engaged but an unmodified eigenvalue limiter is re- 
tained to maintain stability. 

The following defaults will be applied if this hie is not present in the working 
directory: 

If ivisc = 0 in laura_name list .data (see Section 5.4 on page 17) then 
inviscid flow is specified and ivis, jvis, and kvis are set to zero for all blocks 
by default. 

If ivisc = 2 and navier.stokes = .false, in the laura.namelist.data 
hie then the default is a function of the wall boundary condition. If the wall 
boundary on the i min or i max face of block n is a solid surface (see Section 5.3 
on page 14) then ivis(n) = 2, otherwise ivis(n) = 0. In like manner, if the wall 
boundary on the j m i n or j m ax face of block n is a solid surface then jvis(n) = 2, 
otherwise jvis(n) = 0. The default specification for kvis(n) is set to 1 if k mtn or 
k max of block n is a solid surface, otherwise kvis(n) = 0. This default recognizes 
that the standard orientation is for the solid wall on the k m i n surface and the reduced 
eigenvalue limiter is required in this case. 

If ivisc = 2 and navier.stokes = .true, in laura_name list .data then the 
default is ivis = 2, jvis = 2, and kvis = 1 for all the blocks. 

There may be circumstances where the user wishes to override these defaults. If 
a block is away from the stagnation streamline crossing the shock into the boundary 
layer then a more accurate heating on side walls (i and/or j ) will be returned using 
ivis and/or jvis set to 1 without sacrificing stability. In the case of cavities, a 
block may sit over a cavity and not have any solid boundaries itself but has a well 
defined boundary layer streaming into it from an adjacent block. In this case, even 
though the block has no solid boundaries itself, it should engage viscous terms to 
capture the entering shear layer. 

Example: To override the default values and to reset viscous terms in block 

number 3 to ivis = 1, jvis = 0, and kvis = 1, and in block number 5 to ivis =0, 
jvis =1, and kvis =1 the 1 aura_vis.dat a would look like: 
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3 10 1 
5 0 11 


5.12 species_thermo_data 

The species_thermo_data file is the master file for species thermodynamic data. 16 
Each species record consists of the species name, a species properties namelist - 
&species_properties, the number of thermodynamic property curve fit ranges, 
and the curve fit coefficients for each range [15]. No blank line is allowed in this file. 
Elements of the &species_properties namelist are: 

molecule 

A logical flag to determine whether the species is molecule (composed of more 
than one atom); If the species is molecule then molecule = .true., otherwise 
molecule = .false. 

ion 

A logical flag to identify the charged particle, ion = .true, for charged par- 
ticles except for neutrals and electrons. This flag initializes electron-neutral 
energy exchange cross section and sum of the charges. 

charge 

An integer number to determine number of positive charges in particle. If 
ion = .false, then charge = 0. 

elec_impct_ion 

This real number to set the energy for neutrals (i.e. ion=. false.) in electron 
volts (eV) that is required to liberate an electron when the neutral impacts 
with free electron. 

mol_wt 

A real number to set the molecular weight of the particle. This parameter is 
always required. 

siga 

An array of three real numbers, which correspond to curve fit coefficients for 
electron-neutron energy exchange cross section defined as 

a en = a + bT + cT 2 (34) 

where a en is the electron-neutron energy exchange collision cross section in 
m 2 , a, b , and c are the curve fit coefficients, and T is the gas temperature 
[16,17]. The format to define these coefficients is siga=a, b, c. For example, 
siga=7.5e-20, 0, 0. 

16 The species_thermo_data file should only be changed by developers. 
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disoc_ener 


A real number to set dissociation energy of molecule in electron volts (eV). 

alantel 

A real number to set Landau- Teller constant to compute vibrational relaxation 
time for molecule. These are defined in Ref. [18]. This variable is required if 
molecule=.true.. 

cprtO 

A non-dimensional real number that defines translational-rotational heat ca- 
pacity that is normalized by gas constant. This is equal to 

f + 2 

cprtQ = — (35) 

where / is the number of degrees of freedom in translation and rotation. The 
defaults for atoms and diatomic molecules are 2.5 and 3.5, respectively. 


Example: A portion of the species_thermo_data that provides thermodynamic 

properties of carbon is shown below. 


C i 
&species_properties 2 
molecule = .false. 3 
ion = .false. 4 


charge =0 5 
elec_impct_ion = 11.264 6 
siga = 7.5e-20, 5.5e-24, -l.e-28 7 
mol_wt = 12.01070 8 
/ 9 
3 11 


0 . 64950315E+03 
0.19801337E-07 
0 . 85457631E+05 
-0 . 12891365E+06 
0 . 17420927E-06 
0 . 84105978E+05 
0 . 44325280E+09 
0 . 66495953E-06 
0 . 23552734E+07 


-0 . 96490109E+00 
-0.16061440E-10 
0 . 47479243E+01 
0 . 17 195286E+03 
-0.29028178E-10 
0 . 41300474E+01 
-0 . 28860184E+06 
-0.22300788E-10 
-0 . 64051232E+03 


0 . 25046755E+01 
0 . 53144834E-14 
200.000 1000 
0 . 26460444E+01 
0. 16421824E-14 
1000.000 6000 
0 . 77371083E+02 
0.28993887E-15 
6000.000 20000 


-0. 12814480E-04 
0.00000000E+00 
000 

-0.33530690E-03 

0.00000000E+00 

000 

-0.97152819E-02 

0.00000000E+00 

000 


The species name is defined in line 1 . Between lines 2 and 9 species properties are 
defined. These parameters and/or flags state that the carbon molecular weight is 
12.0107, and the species is neither a molecule nor a charged particle, but it can 
liberate an electron when its energy reaches 11.264 eV after it is impacted with a 
free electron. 


Line 10 shows that there are three thermodynamic property curve fits for tem- 
perature ranges of 200 K < T < 1,000 K, 1,000 K < T < 6,000 K, and 6,000 K 


11 

12 

13 

14 

15 

16 

17 

18 
19 
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< T < 20,000 K. Each data range consists of 12 real numbers with a restriction of 
4 real numbers per line. The first 10 real numbers are the thermodynamic curve fit 
coefficients, and the last two real numbers identify the temperature range for the 
given curve fit coefficients. These coefficients are used to calculate the following 
thermodynamic properties 

Cp(T)/R = a\R ~ 1 2 + a 2 T ~ 1 + a 3 + a 4 T + a 5 T 2 + a 6 R 3 + a^T 4 5 6 7 8 ( 36 ) 

h{T)/RT = —a\T 2 + a 2 R 1 ln T + <23 + a 4 — + 05 — + a §— + a 7 ~ — b — ( 37 ) 

rj~ i — 2 JTi 2 rj~\ 3 rj - 4 

s(T)/R = — ai— a 2 T ~ 4 + a^ln T + a 4 T + 05— + ag— + ^7^- + aio ( 38 ) 

where T is the gas temperature, /? is the universal gas constant, c p , h and s are 
the species specific heat, enthalpy and entropy, respectively, and a* are the provided 
curve fit coefficients in species_thermo_data. 

The following corrections will be applied if the gas temperature is out of the 
given range for the given curve fit coefficients: 


Cp(T) = c p (T*) 

h(T) = h(T*) + (T — T*)c p (T*) 
s(T)=s(T*) + ln^c p (T*) 

where 

r~£ j* f Tiower for T T[ ower 

( Ripper fol T > Ripper 

5.13 species_transp_data 

This file 17 contains log-linear curve fit coefficients for species collision cross section 
that are defined based on temperature range of 2,000-4,000 K [19] . This temperature 
range spans boundary-layer temperatures for typical hypersonic entry. The curve fit 
for the given coefficients is poor, however, at temperatures below 1,000 K. A higher 
order curve fit data is available in species_transp_data_0, which supersedes that 
of species _transp_data — see Section 5.14 on the next page. 

1 

2 

3 

4 

5 

6 

7 

8 


Ar Ar 

-14.6017 -14.6502 -14.5501 -14.6028 ! trrl32+kestin et al 
Ar+ Ar+ 

-11.48 -12.08 -11.50 -12.10 
Ar N2 

-14.5995 -14.6475 -14.5480 -14.5981 ! kestin et al 
Ar CO 

-14.5975 -14.6455 -14.5459 -14.5964 ! kestin et al 
17 The species_transp_data file should only be changed by developers. 


(39) 

(40) 

(41) 

(42) 
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5.14 species_transp_data_0 

This file 18 provides collision cross section coefficients for a higher order curve fit 
data [20,21] than those are in the species_transp_data file — see Section 5.13. The 
data in species_transp_data_0 supersedes the data in species_transp_data if the 
file is placed in the working directory hosting the simulation. 

02 N 111 (c) 

-1. 1453028E-03 1 . 2654140E-02 -2 . 2435218E-01 7 . 7201588E+01 

-1.0608832E-03 1 . 1782595E-02 -2 . 1246301E-01 8 . 4561598E+01 

1 . 4943783E-04 -2 . 0389247E-03 1 . 8536165E-02 1 . 0476552E+00 

NO N 111 (c) 

-1 . 5770918E-03 1 . 9578381E-02 -2 . 7873624E-01 9 . 9547944E+01 

-1.47 19259E-03 1 . 8446968E-02 -2 . 6460411E-01 1 . 0911124E+02 

2 . 1014557E-04 -3 . 0420763E-03 2 . 5736958E-02 1 . 0359598E+00 

NO 0 111 (c) 

-1 . 0885815E-03 1 . 1883688E-02 -2 . 1844909E-01 7 . 5512560E+01 

-1.0066279E-03 1 . 1029264E-02 -2 . 0671266E-01 8 . 2644384E+01 

1.4145458E-04 -1 . 9249271E-03 1 . 7785767E-02 1 . 0482162E+00 

END END 110 

0. 0. 0. 0. 

0. 0. 0. 0. 

5.15 surf ace_property_data 

This file is required if there are more than 1 solid surface boundary types (see 
Section 5.3 on page 14) and the surface conditions of these solid surfaces differ from 
those specified in laura_namelist_data (see Section 5.4 on page 17). 

Surface boundary properties of each solid surfaces must be bounded by: 

fesurf ace_properti.es 

/ 

The first instance of the parameters defines properties for surfaces of type - 1 (note 
that properties of type 0 are defined in laura_namelist_data), the second instance 
defines properties for surfaces of type -2, and so on (see Section 5.3 on page 14 for 
different solid surface types.) 

The parameters that can be defined for each solid surface boundary are: 

blowing_model 

See Section 5.4.1 on page 17 for more details and a complete list of options. 
Default: ‘none’ 

18 The species_transp_data_0 file should only be changed by developers. 
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catalysisunodel 

See Section 5.4.10 on page 28 for more details and a complete list of options. 
Default: 'super-catalytic' 

char .density 

Density of the char, kg/rn 3 . Default: 256.29536, which is for the heritage 
AVCOAT. 

CH0NSi_f rac_char 

See Section 5.4.1 on page 17 for more details and a complete list of options. 
Default: CHONSi.f rac.char = 1.0, 0.0, 0.0, 0.0, 0.0 

CH0NSi_f rac_pyrolysis 

See Section 5.4.1 on page 17 for more details and a complete list of options. 
Default: CH0NSi_f rac.pyrolysis = 1.0, 0.0, 0.0, 0.0, 0.0 

emiss.a, emiss_b, emiss_c, emiss.d 

Real numbers to assume emissivity constant coefficients for solid surface bound- 
ary type (see Section 5.3 on page 14) surface emissivity, e, which is calculated 
as 

e = e a + tbT w + e cT%j + (43) 

where T w is the surface temperature. Values for e a , £&, e c , and are defined 
by emiss.a, emiss.b, emiss_c, and emiss.d, respectively. Default: emiss_a= 
0.89, emiss_b=0, emiss_c=0, and emiss_d=0 

h_ablation 

See Section 5.4.1 on page 17 for more details and a complete list of options. 
Default: 0.0 

mdot_pressure 

See Section 5.4.1 on page 17 for more details and a complete list of options. 
Default: 0.0 

surf ace_group_name 

Character descriptor for surfaces with solid surface boundary types — see Sec- 
tion 5.3 on page 14. Any character can be specified to group solid surface 
boundaries. Default: 'undefined surface' 

surf ace.temperature 

Initial wall temperature in K for the solid surface boundary. This variable 
is similar to twall.bc in Section 5.4.10 on page 28. The wall temperature 
stays constant as specified by this parameter if surf ace_temperature_type 
= ‘constant’. Default: value of twall_bc. 

surf ace_temperature_type 

A character identifier for surface temperature model. See Section 5.4.10 on 
page 28 for options and their descriptions. Default: ‘ constant’ 
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t_ablation 


See Section 5.4.1 on page 17 for more details and a complete list of options. 
Default: 0.0 

virgin.density 

Density of virgin material, kg/rn 3 . Default: 544.627742, which is for the 
heritage AVCOAT. 
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6 Output Files 

Laura generates the files listed in Table 2 regardless of inputs specified. A de- 
scription of each output file is presented after a brief discussion of stdout output. 


Table 2: Laura output files. 


Filename 

Content 

Format 

laura.g. fvbnd 

Boundary types in laura.g 

Field View 

laura.nam 

Variables names in laura.q 

Tecplot™ 

laura.q 

Flowfield solution 

Plot3D 

laura_blayer.dat 

Boundary layer edge data 

Tecplot™ 

laura_conv.out 

Time step, CPU time, and residuals 

ASCII 

laura_new.g 

Volume grid 

Plot3D 

laura_new.rst 

Flowfield solution for restart 

Binary 

laura.surf ace.g 

Surface grid 

Plot3D 

laura.surf ace.nam 

Variables names in laura.surf ace.q 

Tecplot™ 

laura.surf ace.q 

Surface solution 

Plot3D 


In addition to files, Laura writes some information on the screen. After initial 
configuration information such as grid block sizes, laura_name list .data specifica- 
tions, and tdata gas physics, the screen output is divided to five distinctive cate- 
gories: block number, task number, L ^ including location and equation, and L 2 of 
mixture continuity, .xyt-momenta, and energy equations. 

The L 2 norm is defined as 


<44) 

where N is the total number of grid cells times the number of conservation equations, 
R is the residual, and p is the mixture density. 

The location of the maximum L ^ residual is shown by four integer numbers after 
the norm itself. The first three integer numbers correspond to i-, j-. and /e-indices of 
the corresponding block, and the last number is the equation number, n_eqn. The 
total number of equations, n_eqn_max, is defined as 

n_eqn_max = n_speci.es + 3 + n.energy (45) 

where n_speci.es and n.energy are number of species and energy equations, respec- 
tively. The equation number corresponding to species is based on the species order 
defined in tdata followed by the x-, y-, and ^-momentum equation numbers. For 
example, consider block 11 for the following sample screen output, 

Reading block sizes from laura.g 

block = 10 task = 10 Linf : 0.63E+02 23 24 86 1 L2eq: 0.31E+00 0.11E-01 0.12E-01 0.12E-01 0.32E-01 0.27E-02 
> block = 11 task = 11 Linf: 0.12E+04 21 2 86 14 L2eq: 0.42E+01 0.10E+01 0.17E+00 0.95E+00 0.88E+00 0.28E-01 < 
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block 

= 

12 

task 

= 

12 

Linf : 

0.31E+03 

11 

7 

87 

5 

L2eq: 

0.20E+01 

0 . 38E+00 

0.22E+00 

0 . 28E+00 

0 . 35E+00 

0.11E-01 

block 

= 

6 

task 

= 

6 

Linf : 

0 . 68E+02 

6 

24 

87 

1 

L2eq: 

0 . 39E+00 

0.17E-01 

0.80E-02 

0.28E-01 

0.37E-01 

0.32E-02 

block 

= 

2 

task 

= 

2 

Linf : 

0 . 66E+02 

3 

24 

87 

1 

L2eq: 

0 . 49E+00 

0.27E-01 

0.88E-02 

0.53E-01 

0.53E-01 

0.38E-02 

block 

= 

4 

task 

= 

4 

Linf : 

0 . 17E+03 

24 

1 

87 

14 

L2eq: 

0 . 89E+00 

0 . 12E+00 

0.47E-01 

0.14E+00 

0.12E+00 

0.75E-02 

block 

= 

9 

task 

= 

9 

Linf : 

0 . 65E+02 

21 

24 

86 

1 

L2eq: 

0 . 43E+00 

0.26E-01 

0.36E-01 

0.47E-01 

0.48E-01 

0.32E-02 

block 

= 

14 

task 

= 

14 

Linf : 

0.78E+02 

2 

24 

86 

1 

L2eq: 

0 . 66E+00 

0.63E-01 

0.27E-01 

0.22E-01 

0.67E-01 

0.36E-02 

block 

= 

7 

task 

= 

7 

Linf : 

0.71E+03 

14 

2 

87 

14 

L2eq: 

0.27E+01 

0.51E+00 

0.11E+00 

0.69E+00 

0 . 55E+00 

0.22E-01 

block 

= 

1 

task 

= 

1 

Linf : 

0 . 22E+03 

24 

4 

87 

14 

L2eq: 

0.11E+01 

0.13E+00 

0.29E-01 

0.14E+00 

0.13E+00 

0.90E-02 

block 

= 

8 

task 

= 

8 

Linf : 

0 . 26E+03 

6 

1 

87 

14 

L2eq: 

0.10E+01 

0 . 14E+00 

0.97E-01 

0.19E+00 

0.16E+00 

0.82E-02 

block 

= 

5 

task 

= 

5 

Linf : 

0.78E+02 

1 

24 

87 

1 

L2eq: 

0.47E+00 

0.22E-01 

0.23E-01 

0.52E-01 

0.50E-01 

0.35E-02 

block 

= 

3 

task 

= 

3 

Linf : 

0 . 45E+03 

24 

16 

87 

14 

L2eq: 

0.21E+01 

0 . 35E+00 

0.80E-01 

0.41E+00 

0 . 34E+00 

0.18E-01 

block 

= 

13 

task 

= 

13 

Linf : 

0.15E+03 

17 

2 

87 

5 

L2eq: 

0 . 98E+00 

0.73E-01 

0.63E-01 

0.65E-01 

0.88E-01 

0.39E-02 

step 

= 

10 

time = 


57.20 

sum (abs (task error)) = 

= 0 

. 28E+02 L2 norm 

= 0.17E-02 





which is for 14- block, 11-species, two-temperature case. The maximum change is 
shown to be from ^-momentum equation n_eqn=14 located on i= 21, j= 2, and k= 86 
of block 1 1 . 

See Section 7.6 on page 58 for a description of the laura_stdout_to_tec utility, 
which can convert this output into a Tecplot™-compatible format. 


6.1 laura.g.f vbnd 

FieldView™ boundary data file used to label boundary condition types contained in 
laura.g file. Tecplot™ automatically detects this file when the laura.g is loaded. 

6.2 laura.nam 

Tecplot™ name file used to label variables contained in laura.q file. 

6.3 laura.q 

Plot 3D function file for post-processing volume solution. Most of the variables 
in this file are non-dimensionalized according to Table 3. The actual number of 
variables in this file depends on the condition specified in the input files. 


Table 3: laura.q variables 


Variables Definition 


Normalized by 


Cn,N 2 ,... 


T 

T 

P 

M u 

P 

c 

e 


Species mass fraction 
Velocity components 
Total energy per unit mass 
Energy mode j per unit mass 
Translational temperature 
Vibrational temperature 
Density 

Molecular weight 
Pressure 
Sonic velocity 

Static energy per unit mass 
Vibrational energy per unit mass 


U n 


' oo 

u 2 

oo 

T 2 

' oo 


ul 


Pc 


PooUl 

C/oo 

u 2 

ul 
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6.4 laura_blayer.dat 

Boundary layer edge, and wall surface properties and shock stand off distance are 
written in this ASCII and Tecplot™ readable file. Below is an example showing 
the header of the file with their orders. Note that the species mass fraction will be 
written in the exact same order as provided by user in tdata file — see Section 5.5 on 
page 35. Also, the vibrational temperature will only be written if two temperature 
is selected in tdata. 

In addition to these properties, angle-of-attack, Mach and Reynolds-nunrber- 
per-grid-unit will be provided in this file as auxiliary parameters. 

TITLE ="BL EDGE PROPERTIES" 

VARIABLES = "xw (m) " 

"yw (m) " 

"zw (m) " 

"rhow (kg/m"3) " 

"pw (Pa) " 

"Tw (K)" 

"Tvw (K)" 

"Hw (J/kg)" 

"muw (Pa.s)" 

"c<sub>N2</sub>w" 

"c<sub>02</sub>w" 

"c<sub>N</ sub>w" 

"c<sub>0</ sub>w" 

"c<sub>NO</sub>w" 

"qw (W/m~2)" 

"tauwx (Pa)" 

"tauwy (Pa) " 

"tauwz (Pa) " 

"rhoe (kg/m~3) " 

"pe (Pa) " 

"Te (K)" 

"Tve (K)" 

"He (J/kg)" 

"ue (m/s) " 

"ve (m/s) " 

"we (m/s) " 

"Me" 

"mue (Pa.s)" 

"c<sub>N2</sub>e" 

"c<sub>02</sub>e" 

" c<sub>N</ sub>e " 

"c<sub>0</ sub>e" 

" c<sub>NO</sub>e " 

"delta (m) " 
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"deltastar (m) " 

"theta (m) " 

"Re-ue (1/m)" 

"CH (kg/m~2 . s) " 

"Stand-off (m) " 

DATASET AUXDATA Common . AngleOf Attack= " -70 . 000 " 
DATASETAUXDATA Common . Ref erenceMachNumber=" 7.130" 
DATASET AUXDATA Common . ReynoldsNumber=" 4037. 

ZONE T="laura_blayer : Block 1" 


6.5 laura.conv.out 


Time steps and residuals history are written in this file. Some of the flow conditions, 
such as angle-of-attack, free stream conditions, Mach number, Reynolds-number- 
per-grid-unit, etc. are also repeated at the beginning of the file. As shown in the 
following sample, step number, clock time, sum of all the residuals in all active tasks, 
and the overall L 2 norm of the residuals are written in this file. The overall L 2 norm 
is defined as: 19 


L2 


T,i(\rhsj\/pi ) 2 

CFL 2 


(46) 


where rh-s is the residual and p is the density. 


Step 
step = 

0 time= 

10 time = 

4.53 

10.87 

sum (abs (task 

error)) = 

0.27E-03 

L2 

norm = 

0.92E-12 

step = 

20 

time = 

17.13 

sum(abs(task 

error)) = 

0.26E-03 

L2 

norm = 

0.88E-12 

step = 

30 

time = 

23.40 

sum (abs (task 

error)) = 

0.25E-03 

L2 

norm = 

0.84E-12 

step = 

40 

time = 

29.65 

sum (abs (task 

error)) = 

0.25E-03 

L2 

norm = 

0.81E-12 


See Section 7.5 on page 58 for a description of the laura_conv_to_tec utility, which 
can convert this file into a Tecplot™-compatible format. 


6.6 laura_new.g 

The Plot 3D grid file includes any changes associated with grid adaptation or grid 
doubling during the run. The name must be changed to laura.g if user wants to 
restart with this new grid on the next run — see Section 5.2 on page 14 for more 
information on the file format. 


6.7 laura_new.rst 

The restart file contains volume and surface data required for restart from end of 
current run. The name must be changed to laura.rst if user wants to restart from 
this new solution file. 

19 The overall L 2 norm definition is different than the L 2 norm defined for each equation 
by Equation 44. 


Users-Manual. tex 


55 


Rev: 43231, November 1, 2009 



6.8 laura_surf ace.g 

Plot3D 3d whole multiblock surface grid file. See Section 5.2 on page 14 for more 
information on the file format. 

6.9 laura.surf ace.nam 

Tecplot™ name file used to label variables contained in laura.surf ace.q file. 

6.10 laura_surf ace.q 

Plot3D function file for post-processing surface solution. The number of variables 
in this file depends on the conditions specified in the laura_namelist_data. Some 
of the surface variables are presented in Table 4. 


Table 4: laura surface. q variables. 


Variables 

Definition 

Unit 

T 

Surface temperature 

K 

P 

Surface pressure 

N / m 2 

7~y> Tz 

Wall shear stresses 

N / m 2 

Qconv 

Convective heat flux 

W/cm 2 

Qrad 

Radiative heat flux 

W / cm 2 

e 

Surface emissivity 

- 

mdot 

Blowing or Suction rate 

kg/ m 2 -s 


Rev: 43231, November 1, 2009 


56 


Users-Manual. tex 



7 Laura Utilities 


Laura has several interactive utilities that automatically generate some of the re- 
quired input files or otherwise aid running and post-processing Laura simulations. 
These utilities are explained here: 


7.1 bounds 

This interactive utility creates laura_bound_data that contains block connectivity 
data. This utility reads the volume grid data from a Plot3D structured grid, 
laura.g. See Section 5.2 on page 14 for more detail on the file format. Here is a 
sample of an interactive session with bounds: 

Enter precision of laura.g : 1 = single, 2 = double 

2 

Do you want all type 9 bounds to default to type 8? 

Enter 1 for yes or 0 for no: (0) 

0 


BLOCK = 1 k = 1 BOUNDARY 

Area = 0 . 26821087E+03 

(Xcntr , Ycntr , Zcntr) = (-0 . 63020306E+02 , 0.91387367E+02,-0.63128987E+03) 
Enter ITYPE: (0) 

0 


The first question is about the precision of laura.g. Note that any grid file 
created by Laura or one of its utilities is created in double-precision. 

The second question is about type 8 and type 9 boundary data. These boundary 
types are given for the block faces that are shared between two blocks. A type 9 
requires an identical orientation of indices across the shared boundary (increasing 
i to increasing i, increasing j to increasing j . and increasing k to increasing k). A 
type 8 accommodates more general connectivity. 


7.2 coarsen 

Use this utility to coarsen the grid file, laura.g, in i-, j-, and/or k-directions. The 
new grid file, laura_new.g, will be double precision regardless of the input grid 
precision. 
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7.3 convert_bouncLdata 

The utility converts old (pre LAURA. 5) STRTf iles/bound_data.strt files to the 
new laura_bound_data format. Usage: convert_bound_data bound_data.strt. 

7.4 convert _1 aura 

This interactive utility converts cases run with prior versions of LAURA. 5 — see 
Appendix A on page 72. This utility generates laura.g, from old.rst, and a new 
restart, laura_new.rst. 

The following file are either required or optional prior to the execution of this 
utility: 

old.2eq 

This file, which may have a different root name, has two-equation turbulent 
model information. This file is optional. 

old.ep+ 

This file, which may have a different root name, has algebraic turbulent model 
information. This file is optional. 

old.qtw 

This file, which may have a different root name, has surface temperature 
information. This file is optional. If this file is provided, free stream density 
and velocity are also needed. Laura uses this information to approximately 
calculate the related parameters needed to be in laura_new.rst. 

old.rst 

This file, which may have a different root name, is the old restart file that 
contains volume and surface data. The utility will ask for the precision of the 
data and the numbers of grid-blocks that are in this file. This file is required. 

laura_bound.dat a 

Use the convert_bound_data utility, or see Appendix A on page 72 to convert 
this file from old format. This file is required. 

7.5 laura_conv_to_tec 

This utility translates the Laura convergence history file, laura_conv.out, de- 
scribed in Section 6.5 on page 55 into a form usable by Tecplot™. 

Usage: laura_conv_to_tec [options] . 

7.6 laura_stdout_to_tec 

To allow task-specific convergence history plotting, this utility translates the Laura 
standard output stream described in Section 6 on page 52 into a convergence history 
file for each task suitable for Tecplot™. 

Usage: laura_stdout_to_tec [options] [file]. 
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7.7 make_assign_tasks 

Given the number and blocks and processes, the make_assign_tasks utility will 
generate a default assign.tasks (point relaxation with k-sweeps). 

Usage: make_assign_tasks n_blocks n_processes. 

7.8 self.start 

This interactive utility generates a single-block structured grid, laura.g, for fam- 
ilies of 2D, axisymmetric and 3D blunt bodies. This utility will also generate 
laura_bound_data. Schematics of several blunt body families are shown in Fig- 
ure 3. Note that for axisymmetric geometries, the symmetry boundaries must be 
on the y-axis. Parameters for defining a probe shape are shown in Figure 4 on the 
next page. 


This interactive utility modifies (shuffles) variables in the Laura restart, laura.rst, 
to enable continuation of simulation if the gas model variables or parameters includ- 
ing number of species, thermal nonequilibrium, radiation, ablation, and turbulence 
model are modified. The users will be prompted to provide necessary information. 



(a) Spherically blunted cone. (b) Custom aerobrake. (c) Asymmetric 3D cone. 


Figure 3: Sample geometries generated by self .start utility. 


7.9 shuf f le.laura 
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Figure 4: Definition of probe shape parameters used by self _start. 
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8 Sample Cases 


8.1 Sphere: 5-species Air, Thermo-chemical Nonequi- 
librium 

A working directory is created in which all requisite files are assembled. The utility 
self_start is used to generate a grid and the initial template for some required 
input files within this working directory. A verbatim transcript of an interactive 
session using self .start follows. 

% self_start 

Select dimensionality: 

1 = Axisymmetric 

2 = Two-dimensional 

3 = Three-dimensional 

1 

Select geometry: 

1 = Conic (cone/wedge, paraboloid, etc.). 

2 = Aerobrake (includes AFE without axis singularity) . 

2 

Select aerobrake type: 

0 = AFE 

1 = hemisphere 

2 = customized aerobrake 

1 

Enter radius (m) { 1.000000}: 

1. 

Select number of cells along symmetry plane. 

30 

Enter number of cells in k direction, prior to any doubling 
16 

At this point the grid file laura.g, the boundary data file laura_bound_data, 
the task assignment file assign.tasks, and the namelist file laura_namelist_data 
are created. The boundary data file does not require any further modification. The 
task assignment file is set for point-implicit relaxation - the standard practice for 
starting any simulation. The namelist file requires editing to define free stream 
conditions and possibly alter default settings. The edited file used for the first run 
of the test case follows. 

&laura_namelist 

velocity_ref = 5000. ! reference velocity, m/s 

density_ref = 0.001 ! reference density, kg/m"3 

tref = 200.0 ! reference temperature, K 

alpha = 0.000 ! pitch angle, degrees 

twall_bc = 500.0 ! initial wall temperature, K 

chem_flag =1 ! 0 chemically frozen, 1 chemical source on 
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therm_flag =1 ! 0 thermally frozen, 1 thermal source on 

irest =0 ! 0 for fresh start, 1 for restart 

ncyc = 2000 ! global steps 

j update = 4 ! steps between update of jacobian 

ntran = 4 ! steps between update of transport properties 

nitfo = 1500 ! number of lst-order relaxation steps 

iterwrt = 200 ! steps between saves of intermediate solution 

rf_inv =3.0 ! inviscid relaxation parameter 

rf_vis =1.0 ! viscous relaxation parameter 

movegrd = 100 ! number of steps between calls to align_shock 

maxmoves = 0 ! maximum number of calls to align_shock 

re_cell =0.1 ! target cell Reynolds number at wall 

fsh =0.6 ! target bow shock position arc length fraction 

kmax_error = 0.005 ! error norm for doubling grid in k-dir. 

kmax_final = 64 ! max number cells in k dir after all doubling 

nexch = 2 ! steps between exchange of info in mpi 

f rac_line_implicit =0.7 ! fraction of line by block tri-dia 

surf ace_temperature_type_0 = 'radiative equilibrium 1 

catalysis_model_0 = 'super-catalytic 1 

emiss_a_0 = 0.89 

ept = 0.010 ! relaxation factor on read eq wall be 

axi_symmetric = .true. 

two_dimensional = .false. 

xmc = 0 . 0000 

ymc = 0.0000 

zme = 0.0000 

grid_conversion_f actor = 1.0000 

sref = 0 . 43633E-01 

cref = 2.0000 

/ 


The first 5 variables of the namelist in this particular template (other variable 
orderings are acceptable) deal with free stream conditions. The user must set these 
values, otherwise the user will get the default values assuming laminar flow of a per- 
fect gas at 5 krn/s and 0.001 kg/rn 3 . In this case both the chem_f lag and therm_f lag 
are reset to 1 to turn on the chemical and thermal source terms. Other template 
values are defined to provide a reasonable compromise between solution robustness 
and convergence rate for a fresh start solution (irest = 0). The template calls for 
2,000 relaxation steps (ncyc) in the initial run with jacobian updates (j update) and 
transport property updates (ntran) requested every four relaxation steps. The first 
1,500 iterations are executed using first-order spatial accuracy (nitfo) — second- 
order accuracy does not contribute significantly to the solution evolution in the 
initial relaxation period. The inviscid and viscous relaxation factors ( rf_inv = 3 
and rf_vis = 1 ) multiply the respective contributions to the Jacobian matrices 
and provide damping of the update. Larger values sometimes improve robustness 
for more energetic flows but are not required in this case. 
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Template values for grid movement are movegrd = 100, maxmoves = 0 (unlimited 
number of moves), re.cell = 0.1, and f sh = 0.6. These values are approximately 
tuned for optimal response in the opening relaxation process from a fresh start 
where the body materializes in a supersonic flow. Allowing the grid to move every 
100 steps provides frequent opportunity to follow the evolution of the shock front 
as it initially is collapsed on the surface and then reflects off the surface into the 
oncoming flow. Setting re.cell to 0.1 provides very tight stretching near the wall. 

If there is a large difference between the wall temperature and the temperature of 
the first cell center off the wall then the upwind algorithm may fail to sense the wall 
- the boundary condition using the Roe’s averaged interface may admit a supersonic 
flow directed toward the wall at the interface. This condition subverts the ability of 
the inviscid, no-slip boundary to properly engage. The tighter near wall resolution 
enables the upwind scheme to sense the wall under all wall temperature conditions 
tested to date. Setting f sh = 0.6 targets the captured bow shock location at 60% of 
the distance between the wall and the inflow boundary, providing adequate margin 
for the shock to reflect outward without striking the inflow boundary prior to the 
next grid update. 

Automated grid doubling in the k direction is controlled by the triggering error 
norm magnitude (kmax.error = 0.005) and the maximum number of cells in the k 
direction (kmax_f inal = 64). Recall that the grid was initialized with only 16 cells 
in the k direction with the self .start utility. The grid will double to 32 cells in the 
k direction (normal to the wall) when the L2 norm first drops below 0.005. It will 
double again when the error norm next falls below this trigger point. The selection 
of a trigger point for more complex problems (three-dimensional, highly energetic) 
may require user experimentation: trigger too high and the grid doubles too early 
causing a lot of extra work; trigger too low and the solution may ring (limit cycle) 
on a coarse grid and never engage a finer grid. 

The template setting therefore updates the boundary condition after completion 
of a forward and backward sweep through the domain. Other template values are not 
discussed here. They are consistent with default values as described in Section 5.4 
and are included for the users convenience. 

Gas property datafiles must now be copied (or linked) from the share/physics ^modules 
installation directory into the working directory prior to executing laura. The files 
include species_thermo_data (thermodynamic curve fit data), species_transp_data 
and species_transp_data_0 (collision cross section data for computing transport 
properties), kinetic.data (chemical kinetic source term data) and tdata (the gas 
model specific to the current simulation). The tdata file is generally the only file 
that requires editing by the user. For the case of 5 species air in thermochemical 
nonequilibrium the file tdata is defined as follows. 

Two Temperature 
N 6.217e-20 
0 7 . 758e-09 

N2 0.737795 
02 0.262205 

NO l.e-09 
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The fresh start run for this example case is now ready to be executed. It is 
assumed that the executable file laura is available in the working directory. This 
case is small enough to be run in interactive mode. For the purposes of discussion, 
it is convenient to capture the output in a file called out_01. 

% laura >& out_01 (for csh) 

% laura > out_01 2>&1 (for sh) 


The user should note several types of information in out_01. The beginning of 
this file contains diagnostic information regarding presence of optional files, free 
stream conditions, and kinetic model diagnostics including warnings regarding ab- 
sence of some allowed third body collision partners. Lines beginning with “step =” 
keep track of the relaxation step number, elapsed wall time, sum of the LI norms 
over all conservation equations, and the L2 norm. Lines beginning with “block 
=” keep track of the LI norm for mixture continuity, x-momentum, y-momentum, 
z-momentum, total energy, and vibrational-electronic energy residuals. The state- 
ment Calling align_shock... appears after every 100 steps (movegrd = 100). The 
statement Saving restart and plot files, followed by intermediate values of 
aerodynamic coefficients appear after every 200 steps (iterwrt = 200). The grid 
doubles automatically after step 472 from 16 to 32 cells in the k direction when the 
error norm first drops below 0.005. 


block = 1 task = 1 err: 0.17E+01 0.30E+00 0.61E-15 0.44E+00 0.38E+00 0.53E-01 

step = 472 time = 24.65 sum(abs(task error)) = 0.29E+01 L2 norm = 0.50E-02 


Increased kmax to 32 

block = 1 task = 1 err: 0.69E+05 0.54E+02 0.31E-14 0.51E+02 0.14E+04 0.40E+03 

step = 476 time = 25.09 sum(abs(task error)) = 0.70E+05 L2 norm = 0.97E+08 


In general, there is a large jump in the error norm following a grid move or grid 
doubling which rapidly diminishes to pre-adjustment levels. A second doubling from 
32 to 64 cells occurs after step 1348. 


block = 1 task = 1 err: 0.15E+01 0.47E+00 0.13E-14 0.34E+00 0.45E+00 0.47E-01 

step = 1348 time = 113.30 sum(abs(task error)) = 0.28E+01 L2 norm = 0.49E-02 


Increased kmax to 64 

block = 1 task = 1 err: 0.47E+05 0.46E+02 0.39E-14 0.41E+02 0.82E+03 0.23E+03 

step = 1352 time = 114.19 sum(abs(task error)) = 0.49E+05 L2 norm = 0.46E+08 


The interim solution for pressure contours after the first 2,000 steps is shown 
in Figure 5a on the facing page. The corresponding surface pressure and heating 
distributions are shown in Figure 5b. The deep blue zone in front of the sphere 
represents undisturbed free stream conditions. The deep red indicated the high 
pressure stagnation region. The captured shock is approximately located at 60% of 
the distance between the spherical surface and the inflow boundary. The surface 
pressures and shock shape will be shown to be nearly converged at this point but 
the heating rates are still far from converged. Converging the boundary layer profile 
is the focus of the remaining relaxation steps. 

Line relaxation is engaged at this point by editing assign.tasks, sweeping 
around the sphere in the i direction and applying line relaxation across the boundary 
layer in the k direction. 
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(a) Pressure contours - interim solution 
after 2000 steps 



z, m 

(b) Surface pressure and heating - interim 
solution after 2000 steps 


Z 



(c) Pressure contours - converged solution 
after 4000 steps. 



(d) Surface pressure and heating - con- 
verged solution after 4000 steps. 



X 


(e) Pressure contours - l.e-12 error norm 
solution after 5600 steps. 



(f) Surface pressure and heating - l.e-12 
error norm solution after 5600 steps. 


Figure 5: Solutions to the sphere test problem: = 5,000 m/s, = 

0.001 kg/m 3 , Tqo = 200 K, 5-species air, thermochemical nonequilibrium, ra- 
diative equilibrium wall, e = 0.89, and super-catalytic wall boundary. 
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The adapted grid and restart files are renamed to start the second run. 

% cp laura_new.g laura.g 
% cp laura_new.rst laura.rst 

Changes or additions to the file laura_namelist_data are indicated below. 

irest = 1 ! 0 for fresh start, 1 for restart 

jupdate = 20 ! steps between update of jacobian 

ntran = 20 ! steps between update of transport properties 

nitfo = 0 ! number of lst-order relaxation steps 

rf_inv =2.0 ! inviscid relaxation parameter 

rf_vis =2.0 ! viscous relaxation parameter 

movegrd = 400 ! number of steps between calls to align_shock 
maxmoves = 3 ! maximum number of calls to align_shock 

f rac_line_implicit = 0.7 

At the conclusion of these next 2,000 steps (4,000) total the L2 norm has dropped 
to 0.74e-06 and the solution is converged. The solution for pressure contours after 
4,000 steps is shown in Figure 5c. The corresponding surface pressure and heating 
distributions are shown in Figure 5d. Significant reduction in heating level and 
smoothing of the stagnation region profile has occurred using the line relaxation 
across the boundary layer during this second set of 2,000 relaxation steps. 

Line relaxation across the entire shock layer (f rac_line_implicit = 1.0) can 
be accommodated in this case if the relaxation factors (rf_inv and rf_vis) are 
increased to 5. The tolerance for convergence of the L2 norm (rrnstol) is set to l.e- 
12. Grid movement is switched off (movegrd = 0). The latest grid and restart files 
are renamed as before to start the third run. The third run reaches the convergence 
criteria in 1,600 additional relaxation steps, a drop of 6 orders of magnitude in 
the L2 norm. The solutions (Figure 5e and Figure 5f) are nearly identical to the 
corresponding figures at 4,000 steps. 


block = 1 

task = 

1 err : 

0 . 14E-04 

0.50E-05 

0. HE- 14 

0.36E-05 

0.40E-05 0.53E-06 

step = 1580 

time = 

335.39 

sum(abs(task error)) 

= 0 . 27E-04 L2 norm 

= 0.10E-11 

block = 1 

task = 

1 err: 

0.14E-04 

0.48E-05 

0. 12E-14 

0.35E-05 

0.39E-05 0.51E-06 

step = 1600 

time = 

339.63 

sum(abs(task error)) 

= 0 . 26E-04 L2 norm 

= 0 . 98E-12 


Aerodynamic Coefficients 

c_x = -0.8854 

c_y = -0.0000 

c_z = 0.6950 

c_l = -0.0000 

c_m = 0.3514 

c_n = 0.0000 


8.2 Coupled radiation procedure 

Starting with the converged non-radiating LAURA solution, the shuf f le_laura 
routine must be applied to the laura.rst file, and the option to convert from un- 
coupled to coupled radiation must be chosen — see Section 5.4.9 on page 28 for more 
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info on radiation flags. The new .rst file created by shuf f le_laura must then 
be renamed to laura.rst. Furthermore, the following lines must be added to the 
laura_namelist_data file: 

radiation = .true, 
nrad = 5000 
frac_rad_new = 0.8 
iinc_rad = 3 
jinc_rad = 3 

The first of these, radiation, must be set to true. The parameter nrad specifies 
the number of flowfield iterations between calls to Hara, f rac_rad_new specifies the 
fraction of the most recent Hara calculation applied to the LAURA solution, and 
iincrrad and jinc_rad, respectively, specify the increment in the treated points 
along the surface and in the spanwise direction. For axisymmetric cases, jinc.rad 
must be 1. For relatively weakly coupled radiation, such as Earth lunar- return 
conditions [22], nrad may be set to 1,500 for the first two calls to Hara, and 3,000- 
5,000 for the subsequent calls. For strongly coupled flows, such as Mars-return entry 
to Earth [11], nrad should be set to 500 for the first two calls to Hara, and 1,000 
for the rest. 

The radiation models applied by the radiation code (HARA) are defined by the 
optional file hara_namelist_data described in Section 5.6 on page 37. If this file 
is not present in the working directory, then the code determines which radiative 
mechanisms to apply based on the species number densities in the flowfield. 


8.3 Unspecified ablation procedure - Coupled 

The recommended procedure for an unspecified ablation computation, meaning the 
ablation rate and wall temperature is computed as part of the flowfield solution 
(instead of being specified by the user), is as follows: 

1: Obtain a non-ablating solution assuming an equilibrium catalytic and radia- 
tive equilibrium wall. Include only species required for a non-ablating solution. 

2: Apply the shuf f le_laura utility to the converged non-ablating solution. Choose 
the ablation option and increase the number of species to the amount required 
to accommodate ablation species. Add the ablation species to tdata. 

3: Modify laura_namelist_data to include the following — see Section 5.4.1 on 
page 17 for more info on ablation parameters: 

surf ace_temperature_type_0 = 'surface energy balance 1 
blowing_model_0 = "equil_char_quasi_steady" 

CH0NSi_f rac_pyrolysis_0 = 0.547, 0.093, 0.341, 0.019, 0.000 
CH0NSi_f rac_char_0 = 0.488, 0.000, 0.273, 0.000, 0.239 

ept = 0.01 
nexch = 2 
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freq_wall = 50 
bprime_flag = 1 
compute_mdot_initial = 1 
ablation_option = 0 
ablation_verbose = .true. 

where CH0NSi_f rac_pyrolysis_0 and CH0NSi_char_pyrolysis_0 should be 
changed to represent the material of interest. Setting bprime_f lag = 1 spec- 
ifies that an approximate film coefficient diffusion model is applied in the 
surface elemental mass balance. This model is robust enough to apply to a 
converged non-ablating flowfield. For this option, freq.wall specifies how 
often the cells at the wall are updated (instead of nexch). A value of 50 seems 
to work for both weakly and strongly ablating cases. In addition, ept rep- 
resents the fraction of the new ablation solution, which includes the ablation 
rate, wall temperature, and wall species. 

4: Run LAURA for roughly 24,000 iterations. During each ablation computa- 
tion, data is printed out for each point on the body, indicated by 1. The 
level of convergence is indicated with mdot residual at the end of ablation 
computation: 

mdot residual = m) 2 (47) 

l 

Usually, mdot residual = l.E-2 or lower indicates that ablation computa- 
tion is adequately converged within 1%. 

5: If the user considers the b-prime approach of sufficient accuracy, then the 
computation is finished. If the user wishes to apply a rigorous diffusion model 
at the surface, consistent with that applied throughout the flowfield, then the 
following modifications should be made to laura_namelist_data: 

freq_wall = 500 
bprime_flag = 0 

Setting bprime_flag = 0 specifies the rigorous diffusion model at the sur- 
face [2]. This model is significantly less robust than the b-prime approach 
(bprime.f lag = 1) , which is why it requires the solution of the b-prime ap- 
proach as an initial condition. With bprime_f lag = 0, the energy equation is 
solved separately from the elemental mass balance and char equilibrium con- 
straints. The number of flowfield iterations between solutions of the energy 
equation is governed by freq_wall, while the other equations are governed 
by nexch. In general, freq_wall should be much greater than nexch. The 
energy equation requires the convective heating, which must be allowed to 
converge to a meaningful value after the wall properties are perturbed. Note 
that with bprime.flag = 0, the ablation calls are significantly quicker than 
for bprime.f lag = 1, and nothing is printed to the screen. 

6: Run Laura until the ablation rate, wall temperature, and convective heating 
reach steady values within one percent. After each solution of the energy 
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equation, an increase in the residual will be seen. Unlike the b-prime approach, 
this value can be reduced within a reasonable number of iterations to le-7, 
or lower. 


8.4 Unspecified ablation procedure - Uncoupled 

The uncoupled ablation analysis (defined in detail in Ref. [2]) differs from the coupled 
analysis in that the influence of ablation on convective heating is treated approxi- 
mately using the blowing correction. In this uncoupled analysis, ablation is never 
introduced into the flowfield, and it consists of simply a post-processing step to the 
non-ablating flowfield. The recommended procedure for an uncoupled ablation anal- 
ysis involves the same first two steps listed in Section 8.3 on page 67 for a coupled 
analysis. After those steps, modify laura_namelist_data to include the following: 

surf ace_temperature_type_0 = 'surface energy balance' 
blowing_model_0 = "equil_char_quasi_steady" 

CH0NSi_f rac_pyrolysis_0 = 0.8822, 0.0283., 0.0866, 0.0029, 0.0 

ept = 1.0 

bprime_flag = 1 

uncoupled_ablation_f lag = 1 

ncyc = 0 ! global steps 

where CH0NSi_frac_pyrolysis_0 and CH0NSi_char_pyrolysis_0, should be changed 
to represent the ablator material of interest. Note that ncyc = 0 is required. The 
final step is to run Laura. This will first call the ablation model. Then it will 
print out an updated laura_surf ace.q file with the computed ablation rate, wall 
temperature, and altered convective heating (from the blowing correction) values. 
Laura will then terminate without executing any flowfield iterations. 
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Appendix A 


Migrating Cases from Prior Versions 

Follow the following steps to start Laura simulations run with Laura prior to 
version 5. 

Step 1. Create a new working directory and change to it 

7. mkdir [Working Directory] 

Step 2. Copy following required files from the old directory: 

7, cp /old/root .rst old.rst 
7, cp /old/root .tsk assign_tasks 

7, cp /old/STRTf iles/bound_data. strt laura_bound_data 

and, depending on the case, these files: 

7, cp /old/root .qtw old.qtw 
7, cp /old/root . 2eq old.2eq 

Step 3. Use the LAURA convert_bound_data utility, or edit laura_bound_data so 
that there are only integer numbers in this file and there is 6 integers 
separated by at least one space per line per computational block. For 
example, the original file may look something like: 

C : CDATADATADATADATADATADATADATADATADATADATADATADATADATADATADATADATA 

C : : : : : : $Name : $ 


data ( itype(i, 
: / 

1), i=l , 6 ) 
1, 1002111, 

2, 

1, 

0 , 

3 / 

data ( itype(i, 2), i=l,6 ) 

: / 1001211, 1, 

2, 

1, 

0 , 

3 / 


C : : CDATADATADATADATADATADATADATADATADATADATADATADATADATADATADATADATA 
Edit the file to look like this: 

1 1002111 2103 

1001211 12103 

Step 4. Run the interactive convert_laura utility and answer all the questions. 

The restart file from LAURA prior to version 5 is single-precision but files 
produced by LAURA 5 and its utilities are all double-precision. This utility 
creates laura.rst and laura.g files. You will have an option either to keep 
the prior-to-version-5 coordinate system orientation or to rotate the grid 
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Figure A6: LAURA coordinate system orientations. 


to the version 5 coordinate system orientation. As of version 5, Laura 
uses + x-axis as the body normal direction while prior to version 5 Laura 
used the +z-axis as the body-normal direction. The two orientations are 
show in Figure A6. 

To run Laura with Laura grid orientation prior to version 5, use the 
following formulas to correct the angle-of-attack a, and the center of the 
moments coordinates: 


Q-new 

— ®-old 90 

(A48) 

(xTYlc) new 

= -(zmc)oi d 

(A49) 

( yTfic) new 

= +{ymc) 0 id 

(A50) 

( zmc) new 

= +(xmc) 0 i d 

(A51) 


Step 5. Copy the example laura_namelist.dat a file to your working directory 

from the [installs.pref ix] /share/laura directory, where install_pref ix 
is the installation prefix specified when Laura was installed. 

% cp [install_pref ix] /share/laura/laura_namelist_data . 

Use this file as a template to define the simulation parameters. Refer 
to Section 5.4 on page 17 for complete list of options. You must change 
irest = 0 to irest = 1 , otherwise the solution will be initialized to free 
stream conditions. 
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Step 6. Edit assign_tasks so that there are only integer numbers in this file. For 
example, the original file may look something like this: 


CASE: 

C : : : : : : $Name : $ 

13 0 10 1 

2 3 0 1 0 2 

blk swp rlx strt stop task 

Edit the file to read like this: 

13 0 10 1 

2 3 0 1 0 2 

Note: you can also use the Laura utility make_assign_tasks to generate 
a new assign_tasks. 

Step 7. If necessary, create the file laura_vis_data (see Section 5.11 on page 45 
for more detail.) 

Step 8. Modify tdata file (see Section 5.5 on page 35) to define the gas model 
condition for your specific simulation. The species order in this file must 
match the species order used in the prior version of Laura. The other 
data files should not be changed. A20 

Step 9. Run Laura as usual — see Section 4 on page 10. 


A20 These files may be changed if different thermodynamic model, curve-fit data, or thermo- 
chemical reaction is needed — see Section 5.5 on page 35 for more detail. 
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Appendix B 


Additional Molecular Band Systems 

This appendix lists molecular band systems available in addition to those listed 
in Section 5.6 on page 37. The band systems listed here are generally weak emitters 
and absorbers, and are therefore not engaged as a default (unlike those listed in 
Section 5.6 on page 37). Therefore, for these band systems to be engaged, the 
following binary flags (on=l/off=0) must be present and set equal to 1 in the 
hara_namelist_data file. 

treat _band_c2_br 

A binary flag activating the C 2 Ballik-Ramsay band system. 
treat _band_c2_da 

A binary flag activating the C 2 Deslandres-d’Azambuja band system. 

treat _band_c2_fh 

A binary flag activating the C 2 Fox-Herzberg band system. 

treat _band_c2_mulliken 

A binary flag activating the C 2 Mulliken band system. 

treat _band_c2_phil ip 

A binary flag activating the C 2 Philips band system, 
treat _band_co3p 

A binary flag activating the CO 3+ band system, 
treat _band_co_angstrom 

A binary flag activating the CO angstrom band system. 

treat _band_co_asundi 

A binary flag activating the CO Asundi band system. 

treat _band_co_triplet 

A binary flag activating the CO triplet band system. 

treat _band_co2 

A binary flag activating the CO 2 band system. 

treat _band_n2_cy 

A binary flag activating the N 2 Carrol-Yoshino band system, 
treat _band_n2_wj 

A binary flag activating the N 2 Worley- Jenkins band system. 
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treat _band_n2_wor ley 

A binary flag activating the N 2 Worley band system, 
treat _band_no_gamma 

A binary flag activating the NO gamma band system. 

treat _band_no_betap 

A binary flag activating the NO beta-prime band system. 

treat _band.no .gamrnap 

A binary flag activating the NO gamma-prime band system. 
treat _band_o2_sr 

A binary flag activating the O 2 Schumann-Runge band system. 

treat. [?] .photo 

A binary flag activating the molecular photo-ionization mechanism [12] for 
[?] specie, where [?] can be o2, n2, and or no. This mechanism is not 
technically a molecular band system. 


Rev: 43231, November 1, 2009 


76 


Users-Manual. tex 



Appendix C 


Trouble Shooting 


This appendix gives some tips if you experience trouble installing or running 
LAURA. This is far from all-inclusive; it is only a modest attempt to capture some of 
the experience shared by customers over the years through the community LAURA- 
users@lists.nasa.gov and private LAURA-support@lists.nasa.gov email lists. 

Of course the first step is to be sure you have the latest version of LAURA — 
releases are announced on the LAURA-news@lists.nasa.gov mailing list. 

C.l Installation 

C.1.1 Unterminated Constant / Line Truncated 

An error during compilation like, 

ifort [...] -DDATADIR=\" [some really long path]\" [. . .] \ 

-c -o datadir_f ile_manager . o datadir_f ile_manager . F90 
In file datadir_f ile_manager . F90 : 30 

search_paths (2) = "[some really long pa 
1 

Error: Unterminated character constant beginning at (1) 

In file datadir_f ile_manager . F90 : 30 
pa 
1 

Warning: Line truncated at (1) 

make [5] : *** [datadir_f ile_manager . o] Error 1 

are due to your installation path violating Fortran’s “free-form” line limit of 132 
characters. Many compilers do not enforce this, but for those that do, compiler 
options are usually available to circumvent this problem — for example, gfortran’s 

-f f ree-line-length-none or g95’s -f free-line-length-huge. 

C.2 Running 

LAURA is designed to recognize when obvious mistakes in input are made: these 
errors will be reported to you via standard output (nominally, the screen). Even 
when the code runs successfully, the user should always check the standard output 
for warning messages. 

The ideas given below are mostly for when the code blows up, and the code’s 
output ’’advice” is not very helpful (NaN detected, segmentation faults, etc.). 

Most of the time, problems running the code can be traced to the use of an 
inadequate grid, block boundary condition errors, or other configuration missteps. 

The grid should look “nice” . There should be “gentle” stretching factors (shoot for 
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less than 1.15), the grid should be as orthogonal as possible — especially in viscous 
regions, and dramatic changes in distance or orientation from one grid line to the 
next should be avoided. Furthermore, the grid minimum spacing (at walls) should 
be appropriate to the problem you are running. In other words, “viscous” grids 
should have reasonably low cell- Reynolds numbers at the wall (re.cell = 0. 1-1.0 
depending on the amount of diffusion); and one should not trying to run inviscid 
flow on a “viscous” grid. 

The first thing to try after a case blows up is to raise the inviscid and viscous 
relaxation factors — see Section 5.4.8 on page 26. Typically, you can start with values 
like rf _inv=4 and rf_vis=3 and end with values like 3 and 1.5, but sometimes you 
need to raise them much higher (e.g., 20 and 10 or 200 and 100) to get a solution 
past a particularly nasty transient. Other things to try include running first-order, 
increasing the frequency of Jacobian, transport, and inter-processor communication 
updates, or successively ramping up the Mach number. 

Running on a coarser version of the grid sometimes helps get the solution going- 
see the coarsen utility described in Section 7.2 on page 57. 

C.2.1 NaNs 

NaNs (Not a Numbers) are a type of IEEE Floating Point Exception. Most compilers 
have options (either during compilation or at runtime via environment variables) to 
trap these exceptions and stop the code in its tracks when one occurs. For example, 
the g95 compiler has environment variables of the form G95_FPU_INVALID that can 
control this while the Intel compiler has the -fpe compilation option. Note: to pin 
point where in the code the exception occurs, you will also want to turn on tracing 
(e.g., -traceback for Intel and -ftrace=full for g95) and symbols (-g) so the 
compiler will give you the precise source-code line number. 

C.2.2 Segmentation Faults 

Segmentation faults usually mean you’ve hit a shell-based memory limit, you’ve run 
out of memory, or you’ve uncovered a coding error. 

To check the first, remove your shell memory limits, i.e., 

csh: limit stacksize unlimited 
bsh: ulimit -s unlimited 
ulimit -d unlimited 
ulimit -m unlimited 

and try re-running. Note: For an MPI job, these have to be in your shell startup 
environment (e.g., .cshrc or .bashrc). 

To check the second, try a smaller case, request more resources (e.g., use less 
cores per node), and/or use the batch option of the top command (-b) to monitor 
memory usage. 

For the third, please see Support, section D on the facing page. 
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Appendix D 


Support 


Because Laura’s primary purpose is to serve NASA missions and it is not offered 
as commercial software, support is available on an “as time allows basis” via two 
channels: the public LAURA-users@lists.nasa.gov community email list, which you 
are encouraged to join via lists.nasa.gov/mailman/listinfo/LAURA-users, and the 
private LAURA-support@lists.nasa.gov email list. 

If your issue is not proprietary or otherwise sensitive, you are strongly encour- 
aged to use and become a member of the LAURA-users@lists.nasa.gov list. 
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